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1.0 INTRODUCTION 


The theory described in Vol. 1 has been Implemented in the MARC nonlinear 
finite element code. Ref. 1. Several modifications to the code are required to 
insure accurate results. The major modification is through the MARC user sub- 
routine HYPELA. Section 2.0 contains a brief description of the MARC code, while 
Section 3.0 discusses the inputs to subroutim HYPELA. Section A.O presents the 
control parameters and data cards required to efficiently Integrate the viscoplastic 
equations presented in Appendix 14 of Vol. 1. 


t 
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2.0 DESCRIPTION OF THE MARC PROGRAM 

The viscoplastic theories presented in Volume 1 have been incorporated into the 
MARC general purpose, nonlinear, finite element program. Ref. 1. This program has 
been developed expressly for nonlinear structural analysis. This computer code 
Involves sophisticated computational algo 4 Ithms and advanced finite element 
formulations, but relies on constitutive models not directly applicable to hot 
section components. 

In this contract, as in the previous contract (Ref. 2), the viscoplastic con- 
stitutive theories were incorporated into the MARC program by means of an initial 
stress technique. All of the material nonlinearity in the constitutive equations 
is incorporated into an initial load vector and treated as a pseudo body force in 
the finite element equilibrium equations. Because the \7iscoplastic constitutive 
theories form a “stiff** system of differential equations, it is necessary to form 
the Incremental constitutive equation appropriate to the finite element load incre- 
ment by meanf of a subincrement technique. 

In the subincrement technique the finite element load increment is split into 
a number of equal subincrements and the viscoplastic constitutive theories are inte- 
grated over the small subincrements to form an accurate representation of the incre- 
mental constitutive equation over the fitiite element load increment. Integration 
over each subincrement was accomplished by means of an explicit Euler forward 
difference method. Provided the subincrements are sufficiently small (so that the 
stability level of the forward difference method is not exceeded), the technique 
has been found to work efficiently and accurately, even for large finite element 
load increments. However, it is difficult for the user to pick efficient subincre- 
mental steps, and there is a considerable incentive to use as few sub increments as 
possible, consistent with the stability of the differential equations comprising 
the constitutive theory. 

In order to understand the flow of information, it is necessary to briefly 
examine the MARC code operation. A summary of the operation of the MARC nonlinear 
finite element, taken from Ref. 2, follows. 

The principle of virtual work may be used to generate the MARC nonlinear 
equilibrium equations governing the Incremental response of the structure to an 
Increment of load. In evaluating the nonlinear structural response of a component, 
the program assumes that the load history is divided into a number of Incrementally 
applied load steps. Each load step is sequentially analyzed as a linear matrix 
problem using an appropriate stiffness matrix and load vector. Although each load 
step uses linear matrix methods to solve the incremental equilibrium equations, the 
incremental equilibrium equations themselves are nonlinear since the load vector 
will depend on the displacement Increment obtained In the solution of the Incre- 
mental equilibrium equations. 
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The principle of virtual work may be written, for applied external point loads 
Pi, or displacements u^, in the form 


Sc] tr, dV:Su]Pj , 


( 1 ) 


where the integral extends over the volume, V, of each finite element and the sum- 
mation sign extends to all the elements in the structure. 

In Eq. (1) the virtual displacement vector 6u^ is related to the virtual 
strain vector 6e^ through the relationship 


8«, = Bij8u,0r8«T=SujB;^, (2) 

Where Bj . is the strain-displacement matrix and the superscript T denotes trans- 
position. Since 6uj^ is an arbitrary virtual displacement vector, Eqs. (1) and 
(2) may be written in the form 






jdV=P| 


(3) 


This relation expresses the equilibrium of the structure when the applied load 
vector is and the stress vector is a^. If an incremental load AP^ is applied 
to the structure and the stress vector changes to + AOj^, the relation expressing 
the equilibrium of the structure at the end of the incremental load application 
may be written as 


l/^B]|(orj+A<Tj)dV:FS+APi 


(A) 


Hence, the relation ecpr ess Ing the equilibrium of the structure during the appHca 
tion of the incremental load vector AP^ is obtained from gqs. (3) and (A) by 
subtraction in the form 


1 JbI Aajdv^AP: 

* J i * 


( 5 ) 
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The MARC code allowf> the user to Implement very general constitutive relation 
ships into the program by means of the user subroutine HYPELA. Within this sub- 
routine the user must specify the values of the elasticity matrix ana “he 
inelastic stress increment vector in the Incremental vector constitutive rela- 
tionship 


AcT|*Djj(A«. -8jO A©) " ■ 


( 6 ) 


The Inelastic stress increment vector is computed in HYPELA from the 
viscoplastic constitutive relationships programmed in the Appendices. 

In Eq. (6) a dnotes the coefficient of thermal expansion and 6j is the 
vector Kronecker delta symbol, 




• ifO< j <3 


'0 if 3 <i <6 


( 7 ) 


For the class of nonlinear visoplastic constitutive relationships under con- 
sideration in this contract, the incremental inelastic stress vector depends 
in a highly nonlinear manner on the Incremental strain vector Acj^. Since Ae^ = 
Auj, the incremental Inelastic stress vector Ac^^ depends in a highly nonlinear 
manner on the nodal displacement vector increment Auj, so that A^^ » A^^ (Au^). 

Substitution of Eq. (6) into (5) produces the incremental equilibrium 
equations for MARC in the form. 


Z Kjj Au.s APj+ARj + I^B^j ACjdv + ZjfB^jSjQA® dv , 


(B) 







where is the elemental elastic stiffness matrix defined by the relation 

dV. (9) 

V ‘ 

The vector AR^ is the residual load correction vector or out-of-equilibrlum 
force vector from the preceding load increment, 




( 10 ) 
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which is added to the current increment in order to restore the structure to 
equilibrium. The nonlinearity in the incremental equilibrium relationship, de- 
fined in Eq. (8), arises because the Inelastic stress Increment vector depends 
nonlinearly on the displacement increment vector Auj . Values of D^j and A^j appro- 
priate to the current Incremental load step are returned to the main program by 
subroutine HYPELA and the Incremental equilibrium relations in Eq. (8) are solved 
by successive iterations. 

The solution of the incremental equilibrium equations in (8) is accomplished 
within the MARC code by the following algorithm. At the start of the increment 
the user subroutine HYPELA is entered to determine the elasticity matrix and 
the incremental Inelastic stress vector A^^. On entry to the subroutine the input 
consists of the strain increment vector ACj^, the temperature Increment A6 , the 
time increment At over which the Incremental external load vector AP^^ is applied 
to the structure, and the values of the stress, strain, temperature and viscoplastic 
state variables at the beginning of the increment. Since the incremental strain 
vector, Aej^ * B^.Auj, can only be accurately determined after the solution to the 
incremental equilibrium relationship in Eq. (8) has yielded the correct incremental 
solution vector Auj , the strain increment vector Ae^ Initially used to generate the 
inelastic stress increment vector A^^ must be estimated. The initial estimate for 
Acj^ is assumed to be the value obtained for Ae^ in the preceding increment. On 
exit from subroutine HYPELA the elasticity matrix D^j and the estimated inelastic 
stress increment vector A^j are passed in<"o the main program. After the values of 
and A^j are obtained for each integration point in the structure, the incre- 
mental equilibrium relationship in Eq. (8) is assembled and solved for the incre- 
mental node displacement vector Auj . The incremental strain vector, Ae^ * BijAuj, 

Is then computed and compared with the initial guess for Aei used to generate 
the inelastic incremental stress vector A^.. If this Incremental strain vector is 
equal, within a user specified tolerance, to the Incremental strain vector used 
to compute ACj in the assembly phase, the solution is assumed to have converged. 
Otherwise, the updated strain Increment vector, obtained from the solution of the 
equilibrium relations in Eq. 8), is passed into subroutine HYPELA, a new vector, 

Ac^ • is computed, and the equilibrium equations resolved to yield an improved 
value of Au^ and Ae^. The process is repeated until the value of the vector Ac^ 
on the assembly phase is equal, within a user specified tolerance, to the value of 
the vector Ae^ on the solution phase. After convergence is achieved, the tempera- 
ture, stress vector, strain vector and vlscoplastlc state variables are updated by 
adding the Incremental values generated during the current Increment to the values 
of these variables at the beginning of the Increment. The program then passes 
on to the next load increment where the process is repeated. A flow chart of the 
Iterative procedure required to Implement vlscoplastlc constitutive theories into 
the MARC program is shown in Fig. 1. 
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3.0 IMPLEMENTATION OF VISCOPLASTIC THEORY IN MARC 
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The Integration of the viscoplastic theory described in Vol. 1 can be introduced 
into the MARC code by means of the user subroutine HYPELA, This routine is called 
at each integration point in each element and supplies the elasticity matrix 
and the Inelastic stress increment vector to the main program. 

The required header cards are: 

SUBROUTINE HYPELA(D,G,E,DE,S,TEMP, DTEMP,NGENS,N,NN,KC,MAT,NDI , 

INSHEAR) 

DIMENSION D(NGENS,NGENS) ,G(NGENS) ,E(NGENS) ,DE(NGENS) ,S(NGENS) 

DIMENSION TEMP(l), DTEMP(l) 

FORTRAN CODING 

RETURN 

END 


where 


D(NGENSjNGENS) is the elasticity matrix D^. defined in this 

subroutine (outout argument) , 


G(NGENS) 

E(NGENS) 

DE(NGENS) 


is the inelastic stress increment vector A;, 
defined in this subroutine (output argument) , 


ft 

is the mechanical strain 

at the beginning of the Incremei t (input argument), 


is the Increment of mechanical strain Ac^ - A^aA0 
(input argument), 


S(NGENS) 


is the stress o. at the beginning of the Increment 
(input argument). 


TEMP(l) is the temperature G at the beginning of the 

increment (input argument) , 

TEMP (2) is the time t at the beginning of the increment 

(input argument) , 


I 
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TEMPO) 

TEMP (4) 
through TEMPO) 


is the cumulative inelastic strain R at the 
beginning of the increment (input argument) , 

are the values of the inelastic strains, c^ 
through cg, at the beginning of the increment 
(input argument) , 


TEMP (10) 
through TEMPOS) 


are the values of the equilibrium stresses 
through at the beginning of the increment 
(input argument). 


TEMP (16) 


is the number of subincrements used to integrate 
the last MARC increment. 


DTEMP(l) 

T)TEMP(2) 

DTE^^P(3) 


is the temperature increment A0 (input argument), 

is the time increment At (input argument) , 

is the increment of cumulative inelastic strain 
AR (output argument) 


DTEMP(4) 
through DTEMPO) 


are the incremental values of Inelastic 
strain Ac^ through Acg (output argument). 


DTEMP(IO) 
through DTEMP(IS) 


are the incremental values of the equilibrium 
stress Afij through Afig (output argument) , 


DTEMP(16) . is output as zero, since the NSPLIT is updated 

automatically. 


NGENS 


N 

NN 

KC 


is the size of the D^. matrix (NGENS ■ 3 for plane 
stress problems, NGENS * 4 for plane strain and 
axisyrametric problems, NGENS ■ 6 for three dimen- 
sional problems) (Input argument), 

is the finite element number (input argument), 

is the integration point number (input argument) , 

is the layer number in shell or beam problems 
(input argument) , 


MAT 


is the material identifier (input argument). 
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NDI is the number of direct stress components (NDl = 2 

for plane stress problems, NDI ■ 3 for plane strain, 
axisymmetric and three dimensional problems) (input 
argument) 

NSHEAR is the number of shear components (NSHEAR « 1 f . 

plane stress, plane strain and axisymmetric problems. 
NSHEAR ~ 3 for three dimensional problems) (input 
argument) 
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4.0 INTEGRATION OF VISCOPLASTIC EQUATIONS 
IN SUBROUTINE HYPELA 


4.1 Self Adaptive Integration 


The values of D 


Ij 


and in the incremental constltutivs relation, 


A<Tj = D jj{ a A0) - A C . . 


( 11 ) 


are obtained by a subincrement method. Incremental values of the variables 
AG, At and (Ae^^ - A^^aAG) for the current finite element load increment are split 
into NSPLIT equal values, and the constitutive equations are Integrated over the 
NSPLIT subincrements to provide accurate values of and ACj^. The subroutine 
HYPELA (Appendix 1) Integrates Walker's vlscoplastic equations and calls subroutine 
HYPCON (Appendix 2) to evaluate the material parameters. HYPCON contains the 
latest estimates for the parameters in the modified Walker's Theory described in 
Vol, 1. Each load increment in a MARC analysis is divided into NSPLIT subincre- 
ments. The integration of the constitutive equations is performed by using forward 
differences with a step size determined by dividing the MARC load increment by 
NSPLIT. Subroutine HYPELA performs the integration in two ways: (1) a fixed step 

size, or (2) a variable step size. In the fixed step size forward difference (KEY 
equal to one at line 32 of HYPELA) , NSPLIT is the same for all MARC load incre- 
ments and sub Increments. 

In the variable step size forward difference (KEY equal to zero in line 32 of 
HYPELA) , NSPLIT is determined by the magnitude of the change in a strain measure 
for every subincrement. The change in the strain measure is defined as 


+ ( 12 ) 

Zfi 

where 

( 13 ) 




•ad 


(14) 
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the quantity Ac Is calculated In line 324 of HYPELA and Is stored as variable 
ERRORO. There are three possible ways to determine NSPLIT. The method depends 
on Che size of ERRORO. If 


ERR0R2 < ERRORO < ERRORl (15) 

Then NSPLIT remains the same for the next sub increment (ERRORl and ERR0R2 
are specified in lines 54 and 55 of HYPELA) . If 

ERRORO < ERR0R2 (16) 

NSPLIT is divided in two for the nixt subincrement and rounded (up) to the 
nearest Integer. If 


ERRORO > ERRORl (17) 

then NSPLIT is doubled and the step is recomputed. The value of NSPLIT at the 
end of the increment is stored in the state variable TEMP(16). The initial value 
of NSPLIT can be specified in an INITIAL STATE BLOCK in the MARC model definition 
cards or in line 31 of HYPELA. The maximum value of NSPLIT is specified by MXSPLT 
(line 57 of HYPELA). If NSPLIT exceeds MXSPLT the message: 

"UNABLE TO REDUCE ERROR IN LESS THAN MXSPLT SUBINCREMENTS" 

is written where the value of MXSPLT is inserted in the WRITE statement. After 
this condition is satisfied the integration is performed using a constant step size. 

4.2 Time Independent Terms and Race of Change of Temperature Term--. 

Two remaining variables must be specified in HYPELA: NONISO (line 64 of 

HYPELA) and INDEP (line 67 of HYPELA). If 

NONISO - 1 (18) 

the change in temperature with respect to time terms will be included. If 

NONISO - 0 (19) 

these terms will not be Included. If 

INDEP « 1 (20) 

the time Independent terms will be included. If 

INDEP • 0 (21) 


these terms will not be included. 
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4.3 Modifications to Other MARC Subroutines 

In order to run viscoplastic xnr*dels in a more efficient manner, the MARC code 
can be modified to run using a constant inverted stiffness matrix. Such a modi- 
fication eliminates the requirement to reassemble and invert the stiffness matrix 
with a subsequent savings in CPU time. This can be accomplished, for example, by: 

(1) changing line 65 of subroutine THRUH (Appendix 3) to 

LOADUQ = 0 (22) 

(2) changing line 116 of subroutine INCDTl (Appendix 4) to 

lASMBL =0, and 

(3) changing line 208 of subroutin INCDTl to 

lASIiBi. = 0 

To insure proper runs only the appropriate lines specifying lASMBL in subroutines 
INCDTl and INCDT2 should be changed. For example, if only BOUNDARY CHANGE cards 
are used in the load incrementation part of the MARC input, then only line 208 

Of INCDTl for specifying lASMBL needs to be changed (along with line 65 of THRUH) • 

Line 753 of the MARC subroutine STEG needs to be changed to that shown in 
Appendix 5 to insure proper running of the MARC code whe i using HYPELA. 

4.4 HYPELA Control Parameters 

The subroutine HYPCON, which calculates the temperature dependent material 
parameters, is called four times by HYPELA at lines 68, 188, 193 and 198. Each of 
these calls evaluates the material parameters at a different temperature. The 
first call determines the elastic constants for whxch the stiffness matrix is 
generated, and will be the same on all increments. The second through the fourth 
calls evaluate the material constants at the median temperature of the subincrement, 
at the beginning of the suoincrement and at the end of the subincrement, respec- 
tively. 

The stiffness matrix is generated with the elastic constants determined by the 
temperature variable SFTEMP (see line 59 of HYPELA) . 

MARC solves the Incremental equilibrium £q. (8) by successive iteratloii. To 
see how the equilibrium equation Iterations are converging one can pick the inte- 
gration point, NPRIN, at line 51 of HYPELA, in element number, NELPR, at line 40 of 
HYPELA, at which the maximum amount of nonlinearity is expected. As subroutine 
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HYPELA is entered on the assembly phase the routine prints out the strain increment 
vector Acj[ and the stress increment vector Ao^ at integration point NPRIN in element 
number NELPR. After the equilibrium equations have been solved for the incremental 
displacement vector subroutine HYPELA is again entered with Ae^ " 

the incremental vectors are printed out on the assembly and solution phase or every 
successive Iteration of the equilibrium equations. In this way the convergence of the 
solution to the incremental equilibrium equations can be followed. If no printout 
of the incremental vectors is required, the variable IPR at line 50 of HYPELA, is 
set equal to zero. 

In order to use subroutine HYPELA, several constants must be defined in the 
subroutine, starting at card number thirty-one (31). These constants are: 


MXSPLT 

S 

maximum number of subincrements allowed, 

NELPR 

S 

element number for printout of incremental stress and strain 
vectors , 

IPR 

s 

1 if stress-strain increment output is required, 

0 if stress-strain increment output is not required 

NPRIN 

:s 

integration point number for printout of incremental 
stress and strain vectors. 

NSPLIT 


number of subincrements per MARC increment. 

SFTEMP 

c 

reference temperature for calculating elastic constants 

ERROR2 

- 

minimum change in strain measure (Eq. 12) allowed before 
a subincrement step size 

ERRORl 

ss 

maximum change in strain measure (Eq. 12) allowed before 
halving subincrement step size 

NONISC 

se 

1 to include rate of change of temperature terms, and equals 
zero if these terms are not included 

INDEP 

m 

1 to include time independent terms, and equals zero if these 
terms are not included. 
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In Appendices 6, 7 and 8 are listings of the data input to simulate the three 
thermomechanical fatigue cycles described in Volume 1. 

4.3 MARC Input Data Set 

In the MARC input data deck two cards are required in the parameter set before 
the END card, lliese are: 


HYPOELAS 

STATE VARS 16 


4 . 6 Temperature Input 


A uniform temperature increment over the structure, together with an appropriate 
time increment, can be specified with the following cards: 

THERIIAL LOADS 

1 , 

5.0, 2.0 

BLANK CARD 

In the above sequence of cards, the first state variable increment of 5.0 refers to 
the uniform temperature increment A© = 5®F over the structure. The second state 
variable increment of 2.0 refers to a time increment of At = 2 seconds. The re- 
mainder of the card and the following BLANK card set the remaining fourteen (14) 
state variable increments to zero. Since the STATE VARS card defines sixteen (16) 
state variables, MARC expects this number as input. However, only the first two 
state variables. A© and At are required as input by HYPELA. The remaining (14) 
state variable Increments are defined within subroutine HYPELA. 

If a nonuniform temperature over the structure is required (the usual case), 
the temperature increments and time increment must be set in user subroutine CREDE. 

This cai> be accomplished with the following header cards: 

SUBROUTINE CREDE (DTDL , M, NSTRES , NEQST , NSTATS ) 

DIMENSION DTDL (NSTATS, NEQST, NSTRES) 

C0MM0N/FAR/DUM,L 

N • (where N •> number of Integration points in element number M) 

DO 2 II • 1, N ' 
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DTDL(1,1,1I) = temperature increment at integration point II 

DTDL(2,1,II) * time Increment (can be made a function of load increment 

number L) 


2 CONTINUE 

RETURN 

END 
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APPENDIX 1. SUBROUTINE HYPELA 


ORIGINAL PAGE 19 
OF POOR QUALITY 


SUBROUTINE HYF£LA(D,G.E,DE,S, TEMP, DTEMP,NGENS,N,NN,KC, MAT, NDI, 
D(NGENS.NGENS).G(NGENS).E(N(XNS).DE(MGENS).S(NGENS} 


SUMSIG(6) ,ET(6) .DSC6) ,ETl 


DIMENSION D(NGENS.NGENS) ,G(NGENS) ,E(NGENS) ,DE(NGENS) ,S(NGENS) 
DIMENSION TmP(l),DTEMP(n 
DIMENSION DSIG(6).D0MEG(6) 

DIMENSION SIG(6).0MEG(6),C(6) 

DIMENSION DC(6).DET(6),0HEGI(6),AB(6) 

DIMENSION SUMSlG(65,ET(6}pSC6),ETl(6) 

DIMENSION DEV(6).OV(6),DAB(6;,c{(6) 

COMMON/AKEV/ravlN 

COMMON/FAR/DUM,INC 

COMMON/ CDC/ ICROUT , NUPTRA , NCPIN.NCYCM. TOLER. XF AC .FLAMB.FRCTOL , MINC , 
t NINST 1 . NINST2 ,NDUM. SCALE .NYD , IKSP , NFCQ , FACS , DISPKE , NCYCLE 
C*****SECOND INVARIANT FdNCTION 

SINV( A ,B , C, D,E ,F)-( A*A+B*B+C*C+2 . * (D*DfE*E+F*F) ) *2 . /3 . 

1234 CONTINUE 

C*****THIS SUBROUTINE RETURNS THE ELASTICITY MATRIX D AND INEIJ^STIC STRESS 
C*****INCREMENT G for the WALKER'S THEORY (B.N.CASSENTI UTRC) 

C*****USING THE DIFFERENTUL FORM OF THE THEORY 
C*****EQUATI0NS are integrated USING A FOWARD DIFFERENCE 
C*****WITH ERROR ESTIMATES FOR REVISING TIME STEPS 
^**** 

C 

C FOR SELF ADAPTIVE INTEGRATION SET KEY-O 
C FOR FIXED STEP INTEGRATION SET KEY-1 AND 

C (1) SET NSPLIT IN TQ1P(16) THROUGH AN INITIAL STATE BLOCK AND/OR 

C (2) SET NSPLIT IN THE IF STATEMENT BELOW 

C 

NSPLIT-TB1P(16) 

IF (NSPLIT.lt. .999)NSPLIT-16 
KEY-1 

IF(N.NE.l) GO TO 7 
IF(NN.NE.t) GO TO 7 
IF( NCYCLE. EQ.O) NWALK-O 
NWALK-NWALK+l 
NQ-NWALK-2*NCYCLE 
NQQ-NCYCLE-1 

C*****DETERMINE IF PLANE STRESS, PLANE STRAIN, AXISYMMETRIC, OR 3-D 
C*****kELTYP-1 for plane STRAIN AND AXISYMMETRIC PROBLEMS 
C*****KELTYP-2 FOR PLANE STRESS PROBLEM 
C*****kelTYP-3 for 3-d PROBLEM 
7 IF(NDI.EQ.3.AND.NSHEAR.EQ.n KELTYP-1 


IF(NDI.EQ.2.AND.NSHEAR.EQ. 11 
IF(NDI.EC.3.AND.NSHEAR.EQ.3J 
C*****SE1 UP CONSTANTS 


KELTYP-1 

XELTYP-2 

KELr/P-3 


IF(INC.EQ.O .AND. ABS(DTEMP(l) ) .LT. l.E-9)DTEMP(l)-l.E-12 

IFaNC.EQ.O .AND. ABS(DTEMP(2) ) .LT. l.E-9)DTEMP(2)-l.E-2 

NELPR-1 

lPR-1 

NPRIN-1 

IF ( (INC+NCYCLE) .NE. 0) NSPL1T-TB1P(16) 

C ERRORl IS MAXIMUM ALLOWABLE STRAIN STEP SIZE 
ERROR l-l.E-4 
ERR0R2-ERR0R1/10. 

C MXSPLT IS THE MAXIMUM VALUE OF NSPLIT ALLOWED 
MXSPLT-65 

XF (KEY .EQ. 1) MXSPLT-NSPLIT 

*§e?m-ISemp - — 

C****NONI^oS**illL*NOf iNCLuiE^OTEMP/D^liERMS ^ 

C****N(X)IS0-1 WILL INTERPOLATE TO INCLUDE DTEMP/DTIME TERMS 
NON IS 0*0 

Ctk***lND£P-0 TIME INDEPENDENT TERMS NOT INCLUDED 
C****INDEP-1 TIME INDEPENDENT TERMS ARE INCLUDED 
IHDEP-l 

CALL HYPCON( INDEP .NONISO , DEGM . 
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70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 

12 

,13 

14 

115 

116 
117 

18 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

135 

136 

lis 

139 

140 

141 

142 


1 EEM . ANUM. AKIM »AK2M .ANIHil. AMM . AM IM ,AN2M .AN3M »AM4M ,AN5M »AN6M ,AN7M , 

2 0ME00M.ARINDM,SIINm.DNl0lM.6N2DTM.DOM]^, 

3 ANM,AliMM,AMUM,ClM.C2M,C3M,C4M,C5M) 

C*****SET ZERO STRAIN INCREMENTS - l.E-8 TO AVOID DIVISION BY ZERO 
DESINV-SINV(DE( 1) ,DE(2) ,DE(3) ,DE(4) ,DE(5) ,DE(6) ) 

DO 1 «J*1 (NGENS 

IF(DESINV.EQ.O.) DE(J)-l.E-8 
1 CONTINUE 

C*****PUX STRESSES AT BEGINNINC OF MARC INCREMENT INTO SIGB ARRAY ACCORDING 
C*****TO ELEMENT TYPE 

GO TO(801,802,803),KELTYP 
801 CONTINUE 
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C*****INITIALIZE STATE VARIABLES ON FIRST ENTRY TO SUBROUTINE. ON SECOND 
C*****AND SUBSEQUENT ENTRIES SKIP INITIALIZATION. 

KEVIN-INC+NCYCLE 
IF(KEVIN.NE.O) GO TO 3 
IF (TBlP(l) .GT. l.E-9) 

IF (TEMP(l) .LT. l.E-9) 

DO 2 J«2,15 
TBIP(J)-l.E-9 

2 CONTINUE 

3 CONTINUE 


SFTEMP-T»IP( 1 ) 
T£MP(l)-SFTEMP 


C*****SET STARTING VALUES OF 
DEG-TEMPO) 

T-TEMP(2) 

R«TcMP(3) 

DO 104 KA*1»6 
J-KA+3 


C(KA, 

lG(KA)-0. 

104 CONTINUE 
KSTEP-0 
ERRORO-0. 

C*****SIART INTEGRATION STEP OVER SUBINCREMENT 
5 SPLIT-NSPLIT 

KSTEP-KSTEP+1 

C WRITE (6,4271) KSTEP.NSPLIT.ERRORO.T 

^I5.5x|'NSPLIt- 


STATE VARIABLES DURING PRESENT MARC INCREMENT 


C4271 FORMAT (IH /KSTEP- 
c iJtme- ^,ireio.3) 

IF (NSPLIT .LE. MXSPLT)G0TO4274 
IF (KEY .EQ. 0) HRITE(6,4273) MXSPLT 
4273 FORMAT 


,I5,5X/ERRORO- %1PE10.3,5X 


4274 


1 ( IH , ' ********************************************************' / 

2 IH /UNABLE TO REDUCE ERROR IN LESS THAN '.15/ SUBINCREMENTS'/ 

3 IH / *************************************4******************' '/ 

KEY-1 


'/) 


IF (ABS(T-(TBIP(2)+DTEMP(2))) .GT. 
SPL1T-DTEMP(2 ) / ( TEMP(2 )+DTEMP( 2)-T: 
IF (ABS(rTIMP(2)+DTEMP(2)-T)/DTEMP( 


C*****seT TEMPERATURE AMD TIME SUl 
6 DDEG-DTEMP( D/SPLIT 
DT-DTEMP12) /SPLIT 
C*****SET ST 


ABS(DTEMP(2)/SPLIT) )G0T06 
>) .LT. l.E-3) G0TO420 


SET STATING VALUES FOR DEVIATO 
PRESS-(SIG(1)'I>SIG(2)+SIG(3))73. 


TORIC STRESSES 
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170 
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216 
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DO 7777 J-1,6 
ALPHA-1. 

IF(J.CT. 3) ALPHA-0. 

7777 

^0’'^ STOAIN INTO ARRAY DET ACCORDING 

C*****EIEMENT TYPE 

111 ?SJ$L^,i»®2.63),KELTYP 
61 CONTINUE 

, > /SPLIT 

2) /SPLIT 

3) /SPLIT 
S*DE(4)/SPLIT 



>-DE(n /SPLIT 
HdE(2) /SPLIT 
— DET(1)-DET(2) 
-0.5*DE(3)/SPLIT 

- 0 . 

0 . 


ETaj-0.5*E(4) 

ET(5)-0. 

ET(6)-0. 

GO TO 71 

62 l»T(r 
DET(2 
DET(3 
DET(4 
DET(5 
DET(6 

ET(l)-E(l) 

ET(2)-E(2) 

ET(3)— ET(1)-ET(2) 

ET(4)-0.5*E(3) 

DET?5)-0. 

DET(6)-0. 

GO TO 71 

63 CONTINUE 

DO 64 J-1,6 
FAC-l. 

IF(J.GT.3)FAO0.5 
DET(J)-FAC*DE(J) /SPLIT 
ET(J)-FAC*E(J) 

64 CONTINUE 
71 CONTINUE 

C*****COrtPUTE TEMPERATURE DEPENDENT MATERIAL CONSTANTS 
DEGM-DEG-K).5*D1£G 
CALL HYPC0N( INDEP . NON ISO , DEGM, 

1 EE.ANU.AKl .AK2 . aAiN.AM.AnI . AIS 2 . AN3 ,AN6 , AN7 , 

2 OM£GO,j&IND,SIlNF.DftlDT,DNiDT,6oMDt, 

3 AN,ALAM,AMU,C1D,C2,C3,C4,C5) 

DEGM-DEG 

CALL HYPC0N(INDEP,N0NIS0,DEGM, 

1 EE0.ANUO.AK1O,AK2O,ANIN0.AMO!aN10.AN2O.AN3O,AN4O,AN5O,AN6O,AN7O, 

2 OMEGOO.AitiNDO SIINPO,DNlDTO,DN2DT6,DOMf)TO, * • » • * 

3 ANTO. ALAMO. AMUO, Cl 0,C20,C30,a0,C50) 

DEGM— DEG4’DD£G 

CALL HYPCON(INDEP.NONISO,DEGM, 

1 EEl,ANUl,AKll,AK2l,ANINl.AMl,ANll,AN2l,AN31,AN41,AN51,AN61,AN71, 

2 OMEG01,AKIND1,SIINF1.DN1DT1.DN2DTI,DOMDT1, 

3 ANl.AlAMl.AMUljXll.ch.C31.C41.C5lj 
c*****sET Initial values of EguiLtBRiuM stress 

DENOM-SINV(C( 1) . C(2) . C(3) , C(4 ) , C(5) , C(6) ) 

DEHOM-DENOM+l.E‘30 

AB( 1 )— OMEGO+2 . *OMEGO* (C( 1 ) *C( 1 )+C(4 ) *C( 4 )+C( 6 ) *C( 6 ) 

1^1 .E^SO/ZDIMOM 

AB( 2 )— OMEGCH-2 . *OMEGO* (C(4 ) *C(4 )+C( 2 ) *C( 2 HC( 5 ) *C( 5 ) 
l+l,E-30)/DEN0M 

AB(3)--^GO+2.*OMEGO*(C(6)*C(6)+C(5)*C(5HC(3)*C(3)+ 

1 1 •E'^BO y /DENOM 

AB(4)-2.*OMEG0*(C(l)*C(4)+C(2)*C(4)+C(5)*C(6)+l.E-30)/ 

IDENOM 

AB(5)-2.*OMEG0*(C(4)*C(6)+C(2)*C(5)+C(3)AC(5)+l.E-30)/ 

IDENOM 

AB(6)-2.*OMEG0*(C(l)*C(6)+C(4)*C(5)+C(3)*C(6)+l.E-30)/ 

IDENOM 
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217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 
229 


231 

232 


234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

m 

278 

279 

280 
281 

iU 

284 

285 

288 

289 

290 


1 J+AB( 2 )+AB( 3 ) 

ALPHA-l. " * 
IF(J.GT.3)ALPHA-0. 

OMEGI ( J)-AB( J )-ALPHA*ABSUM/3 . 
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7124 CONTI 

C*****CCHPUTE DRAG STRESS 

AK-AK1-AK2*EXP(-AN7*R) 

C*****C0HPUTE INELASTIC STRAIN SUBINCREMENTS 
DO 67 J-1,6 

OV( J)- 1 . 5*DEV( J)-0MEG( J) 

67 CONTINUE 

ARG1-SINV(0V( 1 ) ,OV(2) ,0V(3) ,OV(4) ,OV(5) ,OV(6) ) 
ARC1>SQRT(ARG1) 

IF(ABGI.LE.1.E-10)ABG1-1.E-10 
ARG2-(ARGl/AK)**AN 
FAC-(ARC2/ABGl)*DT 
DO 302 J-1 ' 


W — A »W 

DC( J)»FAC* U .5*DEV( J)-C»1EG(J) ) 
302 CONTINUE 

C****ADD TIME INDEPENDENT PART 
IF (INQEF .EQ. 0) G0T068 

DAJ2-1.-AKIND*. ^ 

DWORK-U. 

DO 6801 J-1.6 

6801 


6 666 67*ARG 1* <>2/ (S IINF/AK 1*AK) **2 


Vw w WW L k % Vi 

DWORK-DWORK+SIG( J) *DET( J) *FL0AT( l+J/4) 

IF (DNORK .LT. 0) DWORK-D 

FACT0R-( 1 .-AKIND) *DW0RK/ (S IINF/AK 1*AK) **2/DAJ2 
DO 6803 J-1.6 ^ ^ 

6803 DC(J)-FACT0R*0V(J)+DC(J) 

C*****CC»(PUTE EQUILIBRIUM SIKESS SUBINCREMENTS 
68 OM2«SINV(OMEG( 1) ,aiEG(2) ,0MEG(3) ,0MEG(4),(»1EG(5) ,01EG^6) ) 


(Nt2*'SQRT(0M2) 

DR-SINV(DC( 1) ,DC(2) ,DC(3) ,DC(4) ,DC(5) ,DC(6) ) 

DR-SQRT(DR) 

DG«(AN3+AN4*EXP(-AN5*R) )*DR+AN6*DT*(»!2**(AM-1 . ) 

DO 303 j*1.6 

DOMEG(J)- C ANI+AN2 ) *DC( J)-DG* (OMEG( J)-OMEGI( J)-A1 
IF (NONISO .EQ. 0) GOT0303 . _ 

D0ME^J^-D0MEG^^+(0MEG(J)-0MEGI(J)-AN1*DC(J))*DN2DT*DDEG 


-AN1*C(J)) 




DOM-DGHDT*DDEG 
D(»ll>D0M-2 .*OMEGO*DEENQM 


DENOM-SINV(C( 1 ) , C(2) ,C(3) ,C(4) ,C(5) ,C(6) ) 

DEN0M*DEN0M+1 .E-Bu 

DAB( 1 )— DOlH-2 .*DCM1* (C( 1 )*C( 1 )+C(4) *C(4)+C(6) *C(6) 
141.E-30)7 dENOH 

DAB ( 2 D0l«-2 . *DCM 1 * (C( 4 ) *C(4 HC( 2 ) *C( 2 HC( 5 ) *CC 5 ) 

l-t-l .E-30)/DEN0M 

DAB(3)— DOM+2 .*D011* (C(6) *C(6)+C(5) *C(5)+C(3) *C(3)+ 
^^^4)-2?^0Ml* (C( 1)*C(4)+C(2)*C(4)+C(5)*C(6)+1 .E-30)/ 


IDEN( 

DAB 

IDENi 


DAB^5)-2.*DCMl*(C(4)*C(6HC(2)*C(5)+C(3)*C(5)+l.E-30)/ 


DAB(6)-2.*D(M 1* (C( 1 ) *C(6)+C(4) *C(5)+C(3) *C(6)+1 .E-30)/ 

IDENOM 

DAB( 1 )-^a H2 .*OMEGO* (DC( 1 ) *C( 1 )+DC(4 ) *C(4 )+DC(6) *C(6 ) 

^DAb! 2 2+2 . *0MEG0* (DC(4 ) *C(4 )+DC( 2 ) *C( 2 )+DC( 5) *C( 5 ) 

HUB-SOT/DoiOM 

DAB(3)-DAB(35+2 .*0MEG0* (DC(6)*C(6)+DC(5)*C(5)+DC(3)*C(3)+ 

^DaI(I )-dXb( 4 )+2 . *OMEGO* (DC( 1 ) *C(4 )+DC( 2 ) *C(4 )+DC( 5 ) *C(6 )+ 1 . E-30 ) / 
IDENOK 

DAB« )*DAB( 5 )+2 . *OMEGO* (DC( 4 ) *C( 6 )+DC( 2 ) *C( 5 )+DC( 3) *C( 5 )+l . E-30 ) / 

^DAi[re)-DAB(6)+2.*0MEG0* (DC( 1 )*C(6)+DC(4) *C(5)+DC(3) *C(6 )+l .E-30)/ 
IDLNOM 

DAB ( 1 )-DAB (1 )+2 . *0MEG0* (C( 1 ) *DC( 1 )+C(4 ) *DC( 4 )+C(6 ) *DC( 6 ) 

i+i.b-3o)7den6m 

DAB? 2 )-DAB( 2 )+2 .*0MEG0* (C(4 ) *DC( 4 )+C(2 ) *DC( 2 )+C(5) *DC(5 ) 


20 



.r. - . 




'I 




291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 
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326 
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330 

331 

332 
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334 

335 

336 

337 
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339 

340 

341 
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343 
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l+l.E-30)/DEN0M OF POOR QUALITY 

DAB(3)-DAB(3)+2.*OMEGO*(C(6)*DC(6)+C(5)*DC(5)+C(3)*DC(3)+ 

1 1 *£*3D^/D£M0H 

DAfl(4)-DAB(4)+2.*C»iEG0*{C(l)*DC(4)+C(2)*DC(4)+C(5)*DC(6)+l.E-30)/ 

IDENOM 

DAB(5)-DAB(5)+2.*OMEG0*(C(4)*DC(6)+C(2)*DC(5)+C(3)*DC(5)+l.E-30)/ 

1 DENOH 

nAB(6)»DAB(6)+2.*OMEG0*(C(l)*DC(6)+C(4)*DC(5)+C(3)*DC(6Hl.E-30)/ 

1 DENOM 

303 CONTINUE 

O*****C0MPUTE STOESS SUBINCREMENTS 

IP( KELTYP . EQ . 2) DET( 3 )- ( 2 . * AMU*DC( 3 )- ALAM* (DET( I )+DET( 2) ))/( ALAM+2 . 
1 AAnJ) 

DVOLO-O. 

DVOLl-0. 

DO 809 J-1,6 
ETl(J)-ET(jj+DET(J) 

Cl(J)«C(J)+DC(J) 

IF (J .GT. 3) GOTO809 
DVOLO«DVOLO+ET(J) 

DVOLl-DVOLl+ETl(J) 

809 CONTINUE 

DO 810 J-1,6 
FAC-FLQAT(l-J/4) 

50- ALAMO*DVOLO*FAC+2 .*AMU0* (ET( J)-C( J) ) 

5 1 - ALAM 1 *DVOL 1 *FAC+2 AMU 1 * (E Tl ( J )- C 1 ( J ) ) 

DSIG(J)-S1-S0 

IF (KELTYP ,NE. 2) GOTO810 

DSIGIN-FAO*2 . *AMU*ALAM*DC(3) /(ALAM+2 .*AMU)-2 .*DC(J) *AMU 
DSIG(J)-DSIGIN+( ALAM+2 .*AHJ) *?AC* (DV0Ll-DV0L0)+2 .*AMU*DET( J ) 

810 CONTINUE 

C IF SELF ADAPTIVE INTEGRATION IS USED THEN CALCULATE STRAIN 
C STEP, MAKE APPROPRUTE CHANGES TO NSPLIT AND GO TO 
C CORRESPONDING PROGRAM STEP 
IF (KEY .20. 1)GOT0410 

ERR0R0-SINV(DSIG(1) ,DS1G(2) ,DSIG(3) ,DSIG(4) ,DSIG(5) 

I ,DSIG(6) ) 

ERR0R0-SQRT(ERR0R0)/E&fDR 

IF (ERRORO .LT. ERROR2) NSPLIT-(NSPLIT-1 )/2+l 
IF (ERRORO .LT. ERRORl) G0TO410 
NSPLIT-2*NSPLIT 
G0T05 

C*****UPDATE SUBINCREMENT VARIABLES 
410 DEG-DEG+DDEG 

IF (NSPLIT .GT. MXSPLT) NSPLIT-MXSPLT-1 
R-R-H>R 
DO 113 J-1,6 


»!EG( J )-OMic ( J )+D»!EG ( J ) 
SIG(J)-SIG(J)+DSIG(J) 


Ei(J)-Eii^J) 

SUMSIG(J)-SUMSIG(J)+DSIG(J) 

113 CONTINUE 

PRESS-( S IG( I )+S IG( 2 )+S IG( 3 ) )/ 3 . 

DO 114 J-1,6 
ALPHA-1. 

IF(J.GT.3)ALPHA-0. 
dev1j)-s IG( J)-ALPHA*PRESS 

114 CONTINUE 

C*****END of sub INCREMENT LOOP 
T-T+DT 
G0T05 

C*****PUT elasticity matrix in D and INELASTIC STRESS INCREMENT IN G 
420 GO 10(814,815,816) .KELTYP 
814 CONTI^ 

DO 817 J-l,4 
DO 817 K-l,4 
D(J,K)-0. 

817 CONTINUE 

DO 818 J-1,3 
DO 818 K-1,3 
ALPHA-0. 

IF(J.EQ.K) ALPHA-1. 

D(J^K)-C5M+ALPHA*C3M 
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366 

367 

368 

369 

370 

371 

372 

376 

377 

378 

379 

380 
361 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

il2 

395 

396 

397 

398 

399 

400 

401 

402 
-403 

404 

405 

406 

407 
406 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

ill 

434 

435 

436 


818 CONTI^E 
D(4,4)«C4M 
GO TO 903 


CONTINUE 

D(l,l 

I-C2M 

D(l,2 

I-CIM 

D(2,l 

>-ClM 

Dh,3 

►-0. 

Df3,l 

1-0. 

D(2,2 

►-C2M 

D(2 3 

►-0. 

Dh 2 

>-0. 

DOi3 
GO TO 

I-C4M 

903 
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816 CONTINUE 

DO 819 J-1,6 
DO 819 K«l,6 
D(J,K)-0. 

819 C(^INUE 

DO 820 J-1,3 
DO 820 K-1,3 
ALPHA-0. 

IF(J.EQ.K) ALPKA-1. 
D(J,K)-C5M+ALPHA*C3M 

820 CONTINUE 
D(4,4)-C4M 
D(5,5)-C4M 
Df6.6)-C4M 

903 CONTINUE 

DO 821 J-l.NGENS 
G(J)-SUMSIG(J) 

IF (KELTYP .E' 


^ ^ 2 .AND. J .EQ. 3) G(3)-SUMSIG(4) 

DO 821 K-l.NGgNS 
821 G(J)-G(J)-D(J,K)*DE(K) 

C*****CaMPUTE STRESS AT END OF MARC INCREMENT 
DO 822 J-l.NGENS 
SUM-0. 

DO 823 K-1 .NGENS 
SUM-SUM+DC J , K) *DE(K) 

823 CONTINUE 

DS(j)-SUM-fG(J) 

822 CONTINUE 

C*****PUT STATE VARIABLE INCREMENTS IN TEMP ARRAY FOR NEXT MARC INCREMENT 


DTEMP(3)-R-TEMP(3) 
^ ?(16)-NSPLIT 
DO 923 ibl-1,6 


TEMP( 


J-KA+3 


923 


DTBMP(J)-OMEG(KA)-TEMP(J) 

DTEMPn+6)-C(KA)-TEMP(J+6 

CONTIWE 


) 


iBMAKllH INCREMENT .15) 

E EQ.O) WRITE(6/^3> Noq 
GT.O) WRITE(6.39) NCYCLE 

(55H VALUES OF PARAMETERS DURING SOLUTION OF RECYCLE NUMBER. 


C 

C 


IF(IPR.EQ.O) GO TO 12 
IF(NELPR.NE.N) GO TO 12 
IF (NN.NE.NPRIN) GO TO 12 
WRITE(6,20) INC 
20 FOBMAKllH INCREMENT ^15) 

IF(NQ.EC 
iFr ■ 

23 FOE 

39^FO^T(55H VALUES OF PARAMETERS DURING ASSEMBLY OF RECYCLE NUMBER, 
115) 

WRITE (6, 2 9) 

29 FORMAT(18H STRAIN INCREMENTS) 

WRITE(6,30) (DE(J) ,J-1 , NGENS) 

30 FORMAT(fP6E15.6) 

WRITE(6,3l) 

31 FORMATHSH STRESS INCREMENTS) 

WRITE(6,30) (DS(J).J-l, NGENS) 

12 RETURN 

PTBUG SUBTRACE, UNir(25),SUBCHK,INIT 

AT 1234 

END 
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22 

23 

24 
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SUBROUTINE HYPC0N(INDEP,N0NlS0.DEOl , 

1 EE,ANU,AK1.AK2.AN1N,AM,AN1.AN2,AN3.AN4,AN5,AN6,AN7, 

2 OMEGO , AKIND , SI INF , ON IDT , DN2DT , DOMDT , 

3 AN.ALAM.AMU.C1,C2,C3.C4.C5) 

C THIS SUBROUTINE IS CALLED BY HYPELA TO CALCULATE ALL OF THE 
C TEMPERATURE DEPENDENT MATERIAL CONSTANTS 

DIMENSION TABT(6),EET(6),ANUT(6),AK1T(6),ANINT(6),AMT(6),AN1T(6) 
DIMENSION AN2T(6) ,AN3T(6) ,AN4T(6) ,AN5T(6> ,AN6T(6) ,AN7T(6) ,AK2T(6) 
DIMENSION OMEGOT(^) 

DIMENSION AKINDT(6),SITNFT(6) 

DATA TABT/800. ,1000. ,1200. ,1400. ,1600. ,1800./ 

DATA EET/26.E6,24.E6,24.E6,22.6E6,18.6E6,13.2E6/ 

DATA ANUT/0.32i,0.32d,0.334.0.339’0. 345.6.351/ 

DATA AKlT/5 0931 I .75631. ,9563l. ,ll6696. ,91505. ,59292./ 

DATA AK2T/0. .0. ,0. .0. ,0. ,0./ 

DATA ANINT/.059, .059, .079, .1497, .195. .223/ 

DATA AMT/1.158, 1.158, 1.158, 1.158, 1.158, 1.158/ 

DATA ANlT/0. ,0.,0. ,0.,0. ,0./ 

DATA AN2T/l.E7,l.9E7.1.5E7,2.E7,5.E6,l.E6/ 

DATA AN3T/250. ,320. ,781. 2, 1178. 6, 672. 6, 312. 5/ 

DATA AN4T/0. ,0. ,0. ,0. ,0. ,6./ 

DATA AN5T/0.,0.,0. ,0. 0.,0./ 

DATA AN6T/0. ,0. ,0. ,0. ,8.977E-4,2.733E-3/ 

DATA AN7T/0.,0.,0. ,0.,0. ,0./ 

DATA OMEGOT/0. ,0. ,-2000. ,-2000. ,-1434. ,-1200./ 

DATA AKINDT/0. ,0. ,0. ,0. ,0. ,0./ 

DATA SIINFT/48.E3,48.E3,60.E3,1.E10,1.E10,1.E10/ 

1234 NTP-6 

NTPMl-NTP-1 
TDIF-TABT(2)-TABT( 1 ) 

Ll-DEGM 

L2-TABT(l)-TDIF 

L3-TDIF 

IT-(L1-L2)/L3 

IF(lT.i.T.in>l 

IF( IT. GT . NTPM 1 ) I T^NIPM 1 
FAC-(DEGM-TABT(IT) )/TDIF 
EB-(EET(ITH )-EET( IT) )*FAC+EET( IT) 

ANU-r ANUT( ITU VANUTC IT) )*FAC+ANUT( IT) 

AKl-(AKlTnTfn-AKlTaT) )*FAC+AK1T(IT) 

AK2-(AK2T(lTfl )-^T(IT)>FAC+AK2T(IT) 

ANIN-( AN1NT( ITK> ANINTa T) )*FAC+ANINT( IT) 

AM-(AKr(ITfl)-AMT(IT) )*FAC+AKT(IT) 


ANl- 

AN2- 

AN3- 

AN4- 

AN5- 

AN6» 

AN7- 


ANlT(IT-fl 

AN2T(ITH 

AN3ThTH 

AN4TaTH 

AN5T(ITfl 

AN6T(ITfl 

AN7T(IT+1 


(»1EGO-(OMEGOT( 


I *FAC+AN1T(IT) 
*FAC+AN2T(IT) 
*FAC+AN3T(1T) 
*FAC+AN4T(IT) 
*FAC+AN5T(IT) 
*FAC+AN6T(IT) 
*FAC+AN7T(1T) 

T(IT) )*FAO+OMEGOT(IT) 


-ANIT(IT) 

-AN2T(IT) 

-AN3TfIT) 

-AN4T(IT> 

-AN5TCIT) 

-AN6TfIT) 

-AN7T(IT) 

:Tfl)-0MEG( 

IF (INDFP .EQ. 0) G0T065 

AKIND-TaKINDT( IT f 1 )-AKIl^JDT( IT) )*FAC+AKINDT( IT) 
SIINF-(SIINFT(ITfl)-SIINFTaT))*FAO+SIINFT(IT) 

65 IF (NONISO .EQ. 0) G0T073 

DNlDT-(ANlT(ITf 1 )-ANlT(IT) )/TDIF/ (ANl+l .E-6) 
DN2DT-(AN2T( ITfl )-^2T_(IT) VTDIF/ (AN2*f I .E-6) 
D0MDT-(0MEG6T(ITf 1 )-OMEGOT( IT) )/TDIF/ (OMEGQH .E-6) 
73 CONTINUE 
AN-l./ANIN 

C 1-2 . *AMJ* ALAM/ ( ALAM+2 . * AMO ) 

C2-4 . * AKJ* (ALAM+AMU ) / ( ALAM+2 . * AMU ) 

C3— 2.^AMJ 
C4-AMU 
C5-ALAM 
RETURN 
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DEBUG SUBTRACE , UNIT( 2 5) , SUBCHK » INIT 

AT 1234 

END 
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59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 




lOl 


C 

c * * 


c 

c * * 

c 

c * * 


CALL GMFRD(S,El,DSEG,l,NGENS,l) 

CALL SCLA(B,0. ,NGENS,NGENS,0) 

CALL SCLA(GF.O. .NGENS.I.O) 

USER SUI PLIED HYPOELASTIC CONSTITUTIVE THEORY 


* * * 


CALL HYPELA(£,GF.EELAS,E1,S,DT,0TDL,NGENS.M,NNN,KC,MAT,NDI,NSHEAR) 
LOADUQ*! 

TOTAL ELASTIC STRAINS 

* * * 


CALL GMADD(EELAS,E1,EELAS,NGENS,0 
INCREMENTAL STRESS COMPONENTS 

A * * 
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20A 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 
• 


16 

161 


CONTINUE 

CONTINUE 


c**** boundary change 

IASMBL-1 

READ(5.103)CARD 

CALL INTAC(CARD,ICARD.3,JBAD) 

IF(JBAD.EQ.l)GO TO 110 
NEWBC-ICARD(l) 

IFLAG»ICARD(2) 

IF(IFLAG.EQ.l) WRITE(6,1613) 

1613 FORMAT(63HOFOLLOWING INFORMATION APPLIES TO HARMONIC RESPONSE CALC 
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399 CONTINUE 

U1-MMS4+ISIGXX 

LA2-LMS 1+IGSIG 

IF( I HRESP . EQ . I ) LA2« I GS IG I 

IF(ISHELL.EQ.l)GO TO 617 

IF(JTYPE.EQ.52)GO TO 620 

IF(IHRESP.EQ.l) GO TO 618 

OSEGAA-ABS(OSBG) 

DSEGEA-DSEGE 

OSEGA*DSEG 

IF(ICFEEP.EQ.1.0R.ITHERM.GT.0)CALL TPSMA(VARS(IB1) ,VARS(IECRP1) , 
IVARS(IFCXP),1. ,NGEN1,NGEN1,1) 

IF(IPELA.EQ.l) 

♦CALL GMSUBCVARS(IFCRP),VARS(I<3') .VARS(IFCRP) .NGENl ,1) 
IF(IHERED.EQ.O) GO TO 3l0 
IF(MATV.EQ.O.OR.IVSCFN.NE.2) GOTO 310 
IVS-IVNSER+MATV-1 
NVSER«IOTS(IVS) 

IF(NVSER.EQ.O) GO TO 310 

lVD-IVDSIZ+MATV-1 

NVDSIZ-IKTS(IVD) 

CALL GMADD(VARS(IFCRP),VARS(ICI'V2),VARS(IFCRP),NVDSIZ,l) 

310 CONTINUE 
618 CONTINUE 




■ ^ 


title 

SIZING 

POST 

ALL POINTS 
INPUT TAPE 
HYPOELAS 
STATE VARS 
NO LOADCOR 

restart 

END 

MESH2D 
BLOCKS 
I 4 

DEFINE 
1 2 

boundary 

1 

2 

3 

4 

CONSTRAINT 

2 

I 3 

1 4 

MERGE 
.0005 

generate 

CONNECTIVITY 

1 

COORDINATES 


BOUNDARY CONDITIONS 
0 1 I 

BOUNDARY CONDITIONS 

1 , 

1.3,l,l,l.E-8 
thermal LOADS 
1 , 


KDIX S. DATA INPUT FOR CLOSED SYMMETRIC Uir LOOP 
THER^MECHAN-ICAL loop no. . WITH niNCTIONAL THEORY 
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40 

41 

PROPERTY 

42 

1 

43 

26 • 0E6 

44 

I 4 

45 

POST 

46 

2 16 

47 

1 

48 

11 

49 

PRINT CHOICE 

50 

51 

1 


O.E-5 
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52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 
89 

50 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

!U 

115 

116 

117 

118 

119 

120 
121 


CONTROL 

49 , 15 , 1,0 

.05 

RESTART 
1,1, 0,8,8 
INITIAL STATE 
1.1.1 
1.4,1,4 
800. 

END OPTION 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , .001239 

8.2.1, .001239 

9 . 3 . 1 , .001239 
THERMAL LOADS 
1 , 

3 . 4 , 1.250000 

CONTINUE 
BOUNDARY CHANCE 
3 , 

7 . 1 . 1 , -.000032 

8.2.1, -.000032 

9 . 3 . 1 , -.000032 
THERMAL LOADS 

to. 2, 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000053 

8 . 2 . 1 , -.000053 
9 , 3 ! l " -.000053 
THERMAL LOADS 

{6.8,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000072 

8.2.1, -.000072 

9 . 3 . 1 , -.000072 
THERMAL LOADS 
I, 

23 . 1 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000091 

8.2.1, -.000091 

9 . 3 . 1 , -.000091 
THERMAL LOADS 

19. 1.1.250000 

CONTINUE 
BOUNDARY CHANGE 

?,‘ l , l , -.000108 
8.2 1, -.000108 
9 , 3 , 1 , -.000108 
THERMAL LOADS 

5 4 . 5 . 1.250000 
CONTINUE 
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122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 


BOUNDARY CHANGE 
3, 

7 . 1 . 1 , -.000123 

8 . 2 . 1 , -.000123 

9 . 3 . 1 , -.000123 
THERMAL LOADS 

59.3.1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000136 

8 . 2 . 1 , -.000136 

9 . 3 . 1 , -.000136 
THERMAL LOADS 

43 . 5 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000147 

8 . 2 . 1 , -.000147 

9 . 3 . 1 , -.000147 
THERMAL LOADS 

46 . 9 . 1.250000 

CONTINUE 
BOUNDAT: CHANGE 
3 , 

7 . 1 . 1 , -.000155 

8 . 2 . 1 , -.000155 

9 . 3 . 1 , -.000155 
THERMAL LOADS 

49 . 5 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000160 

8 . 2 . 1 , -.000160 
9 , 3 ! l ! -.000160 
THERMAL LOADS 

51.3.1.250000 

CONTINUE 
BOUNDARY CHANGE 

-.000163 
8 , 2 , 1 , -.000163 

9 . 3 . 1 , -.000163 
THERMAL LOADS 

52 . 2 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 

7 . 1 . 1 , -.000163 

8 . 2 . 1 , -.000163 
9 : 3 , 1 ! -.000163 
TtSRMAL LOADS 
1 , 

52.2,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 
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194 

7,1,1, -.000160 

266 

195 

8,2,1, -.000160 

267 

196 

9,3,1, -.000160 
TmRMAL LOADS 

268 

197 

269 

193 

1, 

51.3,1.250000 

270 

199 

271 

200 


272 

201 

CONTINUE 

m 

202 

BOUNDARY CHANGE 

203 

3. 

275 

204 

7,1,1, -.000155 

276 

205 

8,2,1, -.000155 

277 

206 

9*3!i! -.000155 
THERMAL LOADS 

278 

207 

279 

208 

49.5,1.250000 

280 

209 

281 

210 


282 

211 

CONTINUE 

283 

212 

BOUNDARY CHANGE 

284 

213 

3, 

285 

214 

7,1,1, -.000147 

286 

215 

8,2,1, -.000147 

287 

216 

9,3,1, -.000147 

288 

217 

THERl^L LOADS 

289 

218 

1, 

290 

219 

46.9,1.250000 

291 

220 


292 

221 

CONTINUE 

293 

222 

BOUNDARY CHANGE 

294 

223 

3, 

295 

224 

7,1,1, -.000136 

296 

225 

8,2,1, -.000136 

297 

226 

9,3!i, -,000136 

298 

227 

THERMAL LOADS 

299 

228 

1, 

300 

229 

43.5,1.250000 

301 

230 


302 

231 

CONTINUE 

303 

232 

BOUNDARY CHANGE 

304 

233 

3, 

305 

234 

7ll,l, -.000123 

306 

235 

821 -.000123 

307 

236 

9,3,1, -.000123 

308 

237 

THERMAL LOADS 

309 

238 

1, 

310 

239 

39.3,1.250000 

311 

240 


312 

241 

CONTINUE 

313 

242 

BOUNDARY CHANGE 

314 

243 

3, 

315 

244 

7,1,1, -.000108 

316 

245 

8,2, », -.000108 

317 

246 

9 J 3 I 1 -.000108 

THERMAL LOADS 

318 

247 

319 

248 

1, 

320 

249 

34.5,1.250000 

321 

250 


322 

251 

CONTINUE 

323 

252 

BOUNDARY CHANGE 

324 

253 

3. 

325 

254 

7ll,l, -.000091 

326 

255 

8,2,1, -.000091 

327 

256 

9,3,1, -.000091 
THERMAL LOADS 

328 

257 

329 

258 

1, 

330 

259 

29.1,1.250000 

331 

260 


332 

261 

CONTINUE 

m 

262 

BOUNDARY CHANGE 

263 

3, 

335 

264 

7|l,l, -.000072 

336 

265 

8 2 1 -.0 00072 

337 

338 


-.0 00072 
mRMAL LOADS 

13. 1. 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, -.000053 

8 2^1 -.000053 

9.3.1, -.000053 
t'iermal loads 
I, 

16.8.1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, -.000032 

8.2.1, -.000032 

9!3*1 -.000032 

therMl loads 

10. 2. 1.250000 

CONTINUE 
BOUNDARY CHANCE 
3, 

7.1.1, -.000011 

8.2.1, -.000011 

9.3.1, -.000011 
THE^L LOADS 
1 , 

3.4.1.250000 

CONTINUE 
BOUNDARY CHANGE 

3, 

7.1.1, .000011 

8.2.1, .000011 

9!3!i .000011 

Tten^L LOADS 

1 , 

-3.4,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000032 

8.2.1, .000032 

9.3.1, .000032 
THERMAL LOADS 

-10.2,1.250000 

CONTINUE 
BOUNDARY CHANGE 

3, 

7.1.1, .000053 

8.2.1, .000053 
9’ 3’ I .000053 
THERMAL LOADS 

-16.8,1.250000 

CONTINUE 
BOUNDARY CH/J4GE 
3, 

7.1.1, .000072 
8 2 ll .000072 

9.3.1, .000072 
THEiWAL LOADS 
1 . 
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339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

m 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

¥ii 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 


-23.1,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7.1.1, .000091 

8.2.1, .000091 

9.3.1, .000091 
THERMAL LOADS 

-i9. 1,1.250000 

CONTINUE 
BOUNDARY CHANGE 

?Il,l, .000108 

8.2.1, .000108 

9.3.1, .000108 
THERMAL LOADS 

-54.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000123 

8.2.1, .000123 

9.3.1, .000123 
THEfoiAL LOADS 
1 , 

-59.3,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7;i,l, .000136 

8.2.1, .000136 

9.3.1, .000136 
THEmL LOADS 

-Z3. 5, 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000147 

8.2.1, .000147 

9,3’l, .000147 

THERMAL LOADS 

-46.9,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000155 

8.2.1, .000155 

9.3.1, .000155 
THEmL LOADS 

-49.5,1.250000 


CONTINUE 
BOUNDARY CHANGE 
3, 

7,1,1, .000160 

l±\\ :888II8 

THERMAL LOADS 

1 , 

-il. 3, 1.250000 
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411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 


CONTINUE 
BOUNDARY CHANGE 

711.1, .000163 

8.2.1, .000163 

9.3.1, .000163 
THERMAL LOADS 

-52.2,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000163 

8.2.1, .000163 

9,3,1 .000163 

THEluiAL LOADS 

-52.2,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000160 

8.2.1, .000160 

9.3.1, .000160 
THE^L LOADS 

-51.3,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000155 

8.2.1, .000155 

9’3!i .000155 

THEKIAL LOADS 

-49.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000147 

8.2.1, .000147 

9.3.1, .000147 
THERMAL LOADS 

-46.9,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000136 

8.2.1, .000136 

9.3.1, .000136 
THEK4AL LOADS 

-43.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .0lA.U^ 
8,2 1 .000123 

9.3.1, .000123 
THEmL LOADS 

-59.3,1.250000 

CONTINUE 
BOUNDARY CHANGE 
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483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 


3. 

7.1.1, .000108 

8.2.1, .000108 

9.3.1, .000108 
THERMAL LOADS 

-54.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000091 

8.2.1, .000091 

9.3.1, .000091 
THERMAL LOADS 

-59.1,1,250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000072 

8,2,1 .000072 
9 : 3 ! 1 .000072 

theMl loads 

I, 

-23.1,1.250000 

CONTINUE 
^UNDARY CHANGE 

7, *1,1, .000053 

8,2, ll .000053 

9.3.1, .000053 
THERMAL LOADS 

-16.8,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000032 

8 .2.1, .000032 

9.3.1, .000032 
THERMAL LOADS 

-io. 2, 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3. 

7.1.1, .000011 

8.2.1, .000011 

9.3.1, .000011 
THERMAL IX)A0S 
1 , 

-3.4,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, -.000011 

8 . 2 . 1 , -.000011 
9l3!l’ -.000011 
THERMAL LOADS 
1 , 

3.4,1.250000 

CONTINUE 
BOUNDARY CHANGE 


? 


1 1,1, -.000032 

8,2,1, -.000032 


556 

9,3,1, -.000032 
THEteAL LOADS 

557 

558 

10.2,1.250000 

559 

560 

561 

CONTINUE 

562 

BOUNDARY CHANGE 

563 

3, 

564 

7,1,1, -.000053 

565 

8,2,1, -.000053 

566 

^3’l, -.000053 

567 

THER^ LOADS 

568 

1, 

569 

570 

16.8,1.250000 

571 

CONTINUE 

572 

573 

BOUNDARY CHANGE 
3, 

574 

7|l,l, -.000072 

575 

8,2,1, -.000072 

576 

9,3,1, -.000072 

577 

theMl loads 

578 

1, 

579 

580 

23.1,1.250000 

581 

CONTINUE 

582 

BOUNDARY CHANGE 

583 

3, 

584 

7,1,1, -.000091 

585 

8 2,1. -.000091 

586 

9,3,1, -.0 00091 

587 

THERMAL LOADS 

588 

1, 

589 

590 

29.1,1.250000 

591 

CONTINUE 

592 

BOUNDARY CHANGE 

593 

3, 

594 

7,1,1, -.000108 

595 

8,2,1, -.000108 

596 

9,3,1, -.000108 

597 

THERMAL LOADS 

598 

i4.5. 1.250000 

599 

600 
601 

CONTINUE 

602 

BOUNDARY CHANGE 

603 

3, 

604 

7,1,1, -.000123 

605 

8,2,1, -.000123 

606 

9’3ll -.000123 

607 

THERMAL La\DS 

608 

1, 

609 

610 

39.3,1.250000 

611 

CONTINUE 

612 

BOUNDARY CHANGE 

613 

3, 

614 

7 1,1, -.000136 

615 

8;2,1, -.000136 

616 

9,3,1, -.000136 

617 

THERMAL LOADS 

618 

1 , 

43.5,1.250000 

619 

620 

621 

CONTINUE 

622 

BOUNDARY CHANGE 

623 

3, 

624 

7,1,1, -.000147 

625 

821, -.000147 

626 

9!3,1, -v000147 

627 

THERMAL LOADS 

628 

1. 
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629 

630 

631 

632 

633 

634 

635 

636 

637 

638 

639 

640 

641 

642 

643 

644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

658 

659 

660 
661 
662 

663 

664 

665 

666 

667 

668 

669 

670 

671 

672 

m 

675 

676 

677 

678 

679 

680 
681 
682 

683 

684 

685 

686 

687 

688 

689 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 


46 . 9 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 1 , 1 , -.000155 
8,2 1 . -.000155 

9 . 3 . 1 , -.000155 
THEIMAL LOADS 

49 . 5 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000160 

8 . 2 . 1 , -.000160 
9 ! 3|1 -.000160 
THEMaL LOADS 
1 , 

51 . 3 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000163 

8 . 2 . 1 , -.000163 

9 . 3 . 1 , -.000163 
THERML LOADS 
I, 

52 . 2 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 

3 , 

7 . 1 . 1 , -.000163 

8 . 2 . 1 , -.000163 

9 . 3 . 1 , -.000163 
THERMAL LOADS 

52 . 2 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 

3 , 

7 . 1 . 1 , -.000160 

8 . 2 . 1 , -.000160 

9 . 3 . 1 , -.000160 
THEj^ LOADS 
1 , 

51 . 3 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000155 

8 . 2 . 1 , -.000155 
9 , 3,1 1 -.000155 
THERMAL LOADS 

49 . 5 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 jl , l , -.000147 
8,2 1 , -.000147 

9 . 3 . 1 , -.000147 
THERMAL LOADS 

46 . 9 . 1.250000 


701 

702 

703 

704 

705 

706 

707 

Z 

m 

712 

713 

714 

715 

716 

717 

718 

719 

720 

721 

722 

723 

724 

725 

726 

727 

728 

729 

730 

731 

732 

733 

734 

735 

736 

737 

738 

739 

740 

741 

742 

743 

744 

745 

746 

747 

748 

749 

750 

751 

752 

753 

754 

755 

756 

757 

758 

759 

760 

761 

762 

763 

764 

765 

766 

767 

768 

769 

770 

771 

772 

773 


CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000136 

8 . 2 . 1 , -.000136 

9 . 3 . 1 , -.000136 
THERMAL LOADS 

43 . 5 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 

7 11 . 1 , -.000123 

8 . 2 . 1 , -.000123 

9 . 3 . 1 , -.000123 
TteniAL LOADS 

39 . 3 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000108 

8 . 2 . 1 , -.000108 

9 . 3 . 1 , -.000108 
THEKiAL LOADS 

34 . 5 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000091 

8 . 2 . 1 , -.000091 

9 . 3 . 1 , -.000091 
THERMAL LOADS 

39 . 1 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000072 

8 . 2 . 1 , -.000072 

9 , 3,1 -.000072 

THERMAL LOADS 

33 . 1 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000053 

8 . 2 . 1 , -.000053 

9 . 3 . 1 , -.0 00053 
THERMAL LOADS 

16 . 8 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 

7 . 1 . 1 , -.000032 

8 . 2 . 1 , -.000032 
9 , 3,1 -.000032 
THERMAL LOADS 

10 . 2 . 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3 , 


33 



774 

775 

776 

777 

778 

779 

780 

781 

782 

783 

784 

785 

786 

787 

788 

789 

790 

791 

792 

793 

794 

795 

796 

797 

798 

799 

800 
801 
802 

803 

804 

805 

806 

807 

808 

809 

810 
811 
812 

813 

814 

815 

816 

817 

818 

819 

820 
821 
822 

823 

824 

825 

826 
827 
628 
829 

li? 

832 

111 

835 

836 

837 

838 

839 

840 

841 

842 

843 

844 

845 


7.1.1, -.000011 

8.2.1, -.000011 

9.3.1, -.000011 
THEi^L LOADS 
1 , 

3.4.1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000011 

8.2.1, .000011 

9!3! iI .000011 

THERMAL LOADS 

1 , 

-3.4,1.250000 

CONTINUE 
BOUNDARY CHANGE 

3 

7’, 1,1, .000032 

8,2,1 .000032 

9.3.1, .000032 
THElbiAL LOADS 

-10.2,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000053 

8,2 1, .000053 

9!3|i, .000053 

THE^L LOADS 

-{6.8,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000072 

8 2,1, .000072 

9.3.1, .000072 
THE^L lOADS 

-i3. 1,1.250000 

CONTINUE 
BOUNDARY CHANGE 

3, 

7.1.1, .000091 

8.2.1, .000091 
9*3’ l! .000091 
THERMAL LOADS 

-29.1,1.250000 

CONTINUE 
BOUNDARY CHANGE 

?Il,l, .000108 

6. 2.1, .000108 

9,3.1, .000108 

THERMAL LOADS 

-24.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7 1,1, .000123 

8 2 l| .000123 


846 

847 

848 

849 

850 

851 

852 

853 

854 

855 

856 

857 

858 

859 

860 
861 
862 

863 

864 

865 

866 

867 

868 

869 

870 

871 

872 

873 

874 

875 

876 

877 

878 

879 

880 
881 
882 

883 

884 

885 

886 

887 

888 

889 

890 

891 

892 

893 

894 

895 

896 

897 

898 

899 

900 

901 

902 

903 

904 

905 

906 

907 

908 

909 

910 

911 

912 

111 

915 

916 

917 


9.3.1, .000123 
THERMAL LOADS 
1 , 

-39.3,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000136 

8 2,1 .000136 

9.3.1, .000136 
THERMAL LOADS 

-43.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3. 

7.1.1, .000147 

8;2,1, .000147 

9,3*1 .000147 

THERMAL LOADS 

1 , 

-46.9,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000155 

8.2.1, ,.000155 

9.3.1, .000155 
THERH&L LOADS 

-49.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000160 

8.2.1, .000160 

9!3!i, .000160 

THERMAL LOADS 
1 

-Si. 3, 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7 1,1, .000163 

8.2.1, .000163 

9.3.1, .000163 
THERMAL LOADS 

-32. 2. 1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000163 

8.2.1, .000163 

9l3!l .000163 

THERMAL LOADS 

-32.2,1.250000 

CONTINUE 
BOUNDARY CHANGE 

?Jl,l, ,000160 
8,2,1, .000160 

9.3.1, .000160 
THERMAL LOADS 
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918 

919 

920 

921 

922 

923 

924 

925 

926 

927 

928 

929 

930 

931 

932 

933 

934 

935 

936 

937 

938 

939 

940 

941 

942 

943 

944 

945 

946 

947 

948 

949 

950 

951 

952 

953 

954 

955 

956 

957 

958 

959 

960 
%1 

962 

963 

964 

965 

966 

967 

968 

969 

970 

971 

972 

973 

974 

975 

976 

977 

978 

979 

980 

981 

982 

983 

984 

985 

986 

987 

988 

989 


-51.3,1.250000 

CONTIKUE 
BOUNDARY CHANGE 

3 > 

7.1.1, .000155 

8.2.1, .000155 

9.3.1, .000155 
THE^L LOADS 

-49.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 

^’l,l, .000147 

8.2.1, .000147 

9.3.1, .000147 
THE^L LOADS 

-46.9,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000136 

8.2.1, .000136 

9.3.1, .000136 
THERMAL LOADS 

-43.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 

?Jl,l, .000123 

8.2.1, .000123 

9.3.1, .000123 
THERMAL LOADS 

-59.3.1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000108 

8.2.1, .000108 

9.3.1, .000108 
THERMAL LOADS 

-54.5,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000091 

8.2.1, .000091 

9.3.1, .000091 
THEl^L LOADS 

-19.1,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000072 

8.2.1, .000072 

9.3.1, .000072 
T.WRMAL LOADS 

-13.1,1.250000 


990 

991 

992 

993 

994 

995 

996 

997 

998 

999 
1000 
1001 
1002 

1003 

1004 

1005 

1006 

1007 

1008 

1009 

1010 
1011 
1012 

1013 

1014 

1015 

1016 

1017 

1018 

1019 

1020 
1021 


CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000053 

8.2.1, .000053 

9.3.1, .000053 
THE^L LOADS 

-16.8,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000032 

8.2.1, .000032 

9.3.1, .000032 
THLibiAL LOADS 

-10.2,1.250000 

CONTINUE 
BOUNDARY CHANGE 
3, 

7.1.1, .000011 

8.2.1, .000011 

9,3,1, .000011 

THE^L LOADS 

1 , 

-3.4,1.250000 

CONTINUE 
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APPENDIX 7. DATA INPUT FOR OPEN SYMMETRIC TMF LOOP 


1 

TITLE THERM0ME( 

2 

SIZING 

1500( 

3 

POST 

2 

4 

ALL POINTS 


5 

INPUT TAPE 

1 

6 

HYPOELAS 


7 

STATE VARS 

16 

8 

NO LOADCOR 

1 

9 

RESTART 


10 

END 


11 

MESH2D 


12 

BLOCKS 


13 

1 4 

10 

U 

DEFINE 


15 

1 2 

2 

16 

BOUNDARY 


17 

1 

1 . 

18 

2 

1 . 

19 

3 

0. 

20 

4 


21 

CONSTRAINT 


22 

2 


23 

1 3 

1 

24 

1 4 

2 

25 

MERGE 


26 

.0005 


27 

GENERATE 


28 

CONNECTIVITY 


29 

1 


30 

COORDINATES 


31 


1 

32 

BOUNDARY CONDITIONS 

33 

0 1 

1 

34 

BOUNDARY CONDITIONS 

35 

1, 


36 

1,3,1, 1,1. E- 

8 

37 

THERMAL LOADS 

38 

1 • 


39 

0.,l.E-8 


40 



41 

PROPERTY 


42 

1 


43 

26.0E6 

♦ ' 

44 

1 4 


45 

POST 


46 

2 16 

16 

47 

1 


48 

11 


49 

PRINT CHOICE 


50 


1 


2 WITH FUNCTIONAL THEORY 
2 10 


1 . 

1 . 


I, 




i '■ 


.322 


O.E-5 


36 


lla' 


V', 





51 

1 

104 

7,1,1, -.000126 

52 

CONTROL 

105 

8,2,1, -.000126 

53 

49,15,1,0 

106 

9,3,1, -.000126 

5A 

.05 

107 

THERMAL LOADS 

55 

RESTART 

108 

1, 

56 

1. 1.0,8, 8 

109 

41.1,1.250000 

57 

INITIAL STATE 

no 


58 

M,l 

111 

CONTINUE 

59 

1,4, 1,4 

112 

BOUNDARY CHANGE 

60 

820. 

113 

3, 

61 

END OPTION 

114 

7,1,1, -.000142 

62 

BOUNDARY CHANGE 

115 

8,2,1, -.000142 

63 

3, 

116 

9,3,1, -.000142 

64 

7,1,1, .000150 

117 

THERMAL LOADS 

65 

8,2,1, .000150 

118 

1, 

66 

9,3,1, .000150 

119 

45.0,1.250000 

67 

THERMAL LOADS 

120 


66 

1, 

121 

CONTINUE 

69 

19.0,1 .250000 

122 

BOUNDARY CHANGE 

70 


123 

3, 

71 

CONTINUE 

124 

7,1,1, -.000155 

72 

BOUNDARY CHANGE 

125 

8,2,1, -.000155 

73 

3, 

126 

9,3,1, -.000155 

74 

7,1,1, -.000065 

127 

THERMAL LOADS 

75 

8,2,1, -.000065 

128 

1, 

76 

9,3,1, -.000065 

129 

48.1,1.250000 

77 

THERMAL LOADS 

130 


78 

1, 

131 

CONTINUE 

79 

25.6,1.250000 

132 

BOUNDARY CHANGE 

80 


133 

3, 

81 

CONTINUE 

134 

7,1,1, -.000166 

82 

BOUNDARY CHANGE 

135 

8,2,1, -.000166 

83 

3, 

136 

9,3,1, -.000166 

84 

7,1,1, -.000088 

137 

THERMAL LOADS 

85 

8,2,1, -.000088 

138 

1. 

86 

9,3,1, -.000088 

139 

50.4,1.250000 

87 

THERMAL LOADS 

140 


88 

1, 

141 

CONTINUE 

89 

31.3,1.250000 

142 

BOUNDARY CHANGE 

90 


143 

3, 

91 

CONTINUE 

144 

7,1,1, -.000174 

92 

BOUNDARY CHANGE 

145 

8,2,1, -.000174 

93 

3, 

146 

9,3,1, -.000174 

94 

7,1,1, -.000108 

147 

THERMAL LOADS 

95 

8,2,1, -.000108 

148 

1, 

96 

9,3,1, -.000108 

149 

51.8,1.250000 

97 

THERMAL LOADS 

150 


98 

1, 

151 

CONTINUE 

99 

36.5,1.250000 

152 

BOUNDARY CHANGE 

100 


153 

3, 

101 

CONTINUE 

154 

7,1,1, -.000179 

102 

BOUNDARY CHANGE 

155 

8,2,1, -.000179 

103 

3, 

156 

9,3,1, -.000179 
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157 

THERMAL LOADS 

210 


158 

1, 

211 

CONTINUE 

159 

52.3,1.250000 

212 

BOUNDARY CHA..; £ 

160 


213 

3, 

161 

CONTINUE 

114 

7,1,1, -.000031 

162 

BOUNDARY CHANGE 

215 

8,2.1, -.000031 

163 

3 , 

2Ki 

9,3,1, -.000031 

164 

7,1,1, -.0001 aO 

2' 

THERMAL LOADS 

165 

8,2,1, “-.000180 

218 

1, 

166 

9,3,1, -.000180 

219 

37.5,1.250000 

167 

THERMAL LOADS 

220 


168 

1. 

221 

CONTINUE 

169 

52.0,1.250000 

222 

BOUNDARY CHANGE 

170 


223 

3, 

171 

CONTINITE 

224 

7,1,1, -.000028 

172 

BOUNDARY CHANGE 

225 

8,2,1, -.OO0O28 

173 

3 , 

226 

9,3,1, -.000028 

174 

7,1,1, -.000179 

227 

THERMAL LOADS 

175 

8,2,1, -.000179 

228 

1, 

176 

9,3,1, -.000179 

229 

32.4,1.250000 

177 

THERMAL LOADS 

230 


178 

1, 

231 

CONTINUE 

179 

50.7,1.250000 

232 

BOUNDARY CH/J^GE 

180 


233 

3, 

181 

CONTINUE 

234 

7,1,1, -.000024 

182 

BOUNDARY CHANGE 

2:5 

8,2,1, -.000024 

183 

3 1 

236 

9,3,1, -.000024 

184 

7,1,1, -.000175 

237 

THERMAL LOADS 

185 

8,2,1, -.000175 

238 

1, 

186 

9,3,1: -.000175 

235 

26.8, 1.250000 

187 

THERMAL LOADS 

240 


188 

1. 

241 

C'^NTINUE 

189 

48.0,1.250000 

242 

BuUNDARY CHANGE 

190 


243 

3, 

191 

CONTINUE 

244 

7,1,., -.000020 

192 

BOUNDARY CHANGE 

245 

8, ',1, -.000020 

193 

3 , 

246 

9,3,1, -.000020 

194 

7,1,1, -.000243 

247 

THERMAL lOADS 

195 

8,2,1, -.000243 

2<>S 

1, 

196 

9,3,1, -.000243 

249 

20.^,1.250000 

197 

THERMAL LOADS 

250 


198 

1, 

251 

CONTINUE 

199 

45.7,1.250000 

252 

BOUNDARY CHANGE 

200 


253 

3, 

201 

CONTINUE 

254 

7,1,1, -.000015 

202 

BOUNDARY CHANGE 

255 

8,2,1, -.000015 

203 

•» 

1 

256 

9,3,1, -.000015 

204 

7,1,1, -.000034 

257 

THERMAL LOADS 

205 

8,2,1, -.000034 

258 

i y 

206 

9,3,1, -.000034 

259 

14.2,1.250000 

207 

THER.MAL LOADS 

260 


208 

1 , 

261 

CONTINUE 

209 

41.9,1.250000 

262 

BOUNDARY CHANGE 
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263 

3, 

316 

9,3,1, .000086 

264 

7,1,1, -.000011 

317 

THERMAL LOADS 

265 

8,2,1, -.000011 

318 

1, 

266 

9,3,1, -.000011 

319 

-25.6,1.250000 

267 

THERMAL LOADS 

320 


268 

1, 

321 

CONTINUE 

269 

7.5,1.250000 

32'' 

BOUNDARY CHANGE 

270 


323 

3. 

271 

CONTINUE 

324 

7,1,1, .000114 

272 

BOUNDARY CHANGE 

325 

8,2,1, .000114 

273 

3, 

326 

9,3,1, .000114 

274 

7,1,1, .000108 

327 

THERMAL LOADS 

275 

8,2,1, .000108 

328 

\ , 

276 

9,3,1, .000108 

329 

-31.3,1.250000 

277 

THERMAL LOADS 

330 


278 

1 

- y 

331 

CONTINUE 

279 

.7,1.250000 

332 

BOUNDARY CHANGE 

280 


333 

3, 

281 

CONTINUE 

334 

7,1,1, .000139 

282 

BOUNDARY CHANGE 

335 

8,2,1, .C00139 

283 

3, 

336 

9,3,1, .000139 

284 

7,1,1, -.00:063 

337 

THERMAL LOADS 

285 

8,2,1, -.000063 

338 

1, 

286 

9,3,1, -.000063 

339 

-36.5,1.250000 

287 

THERMAL LOADS 

340 


288 

1, 

341 

CONTINUE 

289 

-6.1,1.250000 

342 

BOUNDARY CHANGE 

290 


343 

3. 

291 

CONTINUE 

344 

7,1,1, .000162 

292 

BOUNDARY CHANGE 

345 

8,2,1, .000162 

293 

3, 

346 

9,3,1, .000162 

294 

7,1,1, .000569 

347 

THERMAL LOADS 

295 

8,2,1, .000569 

348 

1, 

296 

9,3,1, .000569 

349 

-41.1,1.250000 

297 

THERMAL LOADS 

350 


298 

1, 

351 

CONTINUE 

299 

-12.9,1.250000 

352 

BOUNDARY CHANGE 

300 


353 

3, 

301 

CONTINUE 

354 

7,1,1, .000183 

302 

BOUNDARY CHANGE 

355 

8,2,1, .000183 

303 

3, 

356 

9,3,1, .000183 

304 

7,1,1, .001229 

357 

THERMAL J.0ADS 

305 

8,2,1, .001229 

358 

1, 

306 

9,3,1, .001229 

359 

-45.0,1.250000 

307 

THERMAL LOADS 

360 


308 

1, 

361 

CONTINUE 

309 

-19.4,1.250000 

362 

BOUNDARY CHANGE 

310 


363 

3, 

311 

CONTINUE 

364 

7,1,1, .000200 

312 

BOUNDARY CHANGE 

365 

8,2,1, .000200 

313 

3, 

366 

9,3,1, .000200 

314 

7,1,1, .000086 

367 

THERMAL LOADS 

315 

8,2,1, .000086 

368 

1, 
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369 

-48.1,1.250000 

422 

BOUNDARY CHANGE 

370 


423 

3, 

371 

CONTINUE 

424 

7,1,1, .000047 

372 

BOUNDARY CHANGE 

425 

8,2,1, .000047 

373 

3, 

426 

9,3,1, .000047 

374 

7,i,l, .000214 

427 

THERMAL LOADS 

375 

8,2,1, .000214 

428 

1, 

376 

9,3,1, .000214 

429 

-48.6,1.250000 

377 

THERMAL LOADS 

430 


378 

J , 

431 

CONTINUE 

379 

-50.4,1.250000 

432 

BOUNDARY CHANGE 

380 


433 

3, 

381 

CONTINUE 

434 

7,1,1, .000045 

382 

BOUNDARY CHANGE 

435 

8,2,1, .000045 

383 

3, 

436 

9,3,1, .000045 

384 

7,1,1, .000224 

437 

THERMAL LOADS 

335 

8,2,1, .000224 

438 

1. 

386 

9,3,1, .000224 

439 

-45.7,1.250000 

387 

THERMAL LOADS 

440 


388 

1, 

441 

CONTINUE 

389 

-51.8,1.250000 

442 

BOUNDARY CHANGE 

390 


443 

3, 

391 

CONTINUE 

444 

7,1,1, .000043 

392 

BOUNDARY CHANGE 

445 

8,2,1, .000043 

393 

3, 

446 

9,3,1, .000043 

394 

7,1,1, .000230 

447 

THERMAL LOADS 

395 

8,2,1, .000230 

448 


396 

9,3,1, .000230 

449 

-41.9,1.250000 

397 

THERMAL LOADS 

450 


398 

1, 

451 

CONTINUE 

399 

-52.3,1.250000 

452 

BOUNDARY CHANGE 

400 


453 

3, 

401 

CONTINUE 

454 

7,1,1, .000039 

402 

BOUNDARY CHANGE 

455 

8,2,1, .000039 

403 

3, 

456 

9,3,1, .000039 

404 

7,1,1, .000247 

457 

THERMAL LOADS 

405 

8,2,1, .000247 

458 

1. 

406 

9,3,1, .000247 

459 

-37.5,1.250000 

407 

THERMAL LOADS 

460 


408 

1, 

461 

CONTINUE 

409 

-52.0,1.250000 

462 

BOUNDARY CH/JIGE 

410 


463 

3, 

411 

CONTINUE 

464 

7,1,1, .000035 

412 

BOUNDARY CHANGE 

465 

8,2,1, .000035 

413 

3, 

466 

9,3,1, .000035 

414 

7,1,1, .000048 

467 

THERMAL LOADS 

415 

8,2,1, .000048 

468 

1, 

416 

9,3,1, .000048 

469 

-32.4,1.250000 

417 

THERMAL LOADS 

470 


418 

1, 

471 

CONTINUE 

419 

-50.7,1.250000 

472 

BOUNDARY CHANGE 

420 


473 

3, 

421 

CONTINUE 

474 

7,1,1, .000030 



475 

8,2,1, .000030 

531 

CONTINUE 

476 

9,3,1, .000030 

532 

BOUNDARY CHANGE 

477 

THERMAL LOADS 

533 

3, 

478 

1, 

534 

7,1,1, -.000569 

479 

-26.8,1.250000 

535 

8,2,1, -.000569 

480 


536 

9,3,1, -.000569 

481 

CONTINUE 

537 

THERMAL LOADS 

482 

BOUNDARY CHANGE 

538 

u 

483 

3, 

539 

12.9,1.250000 

484 

7,1,1, .000025 

540 


485 

8,2,1, .000025 

541 

CONTINUE 

486 

9,3,1, .000025 

542 

BOUNDARY CHANGE 

487 

THERMAL LOADS 

543 

3, 

488 

1. 

544 

7,1,1, -.001229 

489 

-20.7,1.250000 

545 

6,2,1, -.001229 

490 


546 

9,3,1, -.001229 

491 

CONTINUE 

547 

THERMAL LOADS 

492 

BOUNDARY CHANGE 

548 

1, 

493 

3, 

549 

19.4,1.250000 

494 

7,1,1, .000019 

550 


495 

8,2,1, .000019 

551 

CONTINUE 

496 

9,3,1, .000019 

552 

BOUNDARY CHANGE 

497 

THERMAL LOADS 

553 

3, 

498 

1. 

554 

7,1,1, -.000067 

499 

-14.2,1.250000 

555 

8,2,1, -.000067 

500 


556 

9,3,1, -.000067 

501 

CONTINUE 

557 

THERMAL LOADS 

502 

BOUNDARY CHANGE 

558 

1, 

503 

3, 

559 

25.6,1.250000 

504 

7,1,1, .000013 

560 


505 

8,2,1, .000013 

561 

CONTINUE 

506 

9,3,1, .000013 

562 

BOUNDARY CHANGE 

507 

THERMAL LOADS 

563 

3, 

508 

1, 

564 

7,1,1, -.000088 

509 

-7.5,1.250000 

565 

8,2,1, -.000088 

510 


566 

9,3,1, -.000088 

511 

CONTINUE 

567 

THERMAL LOADS 

512 

BOUNDARY CHANGE 

568 

1, 

513 

3, 

569 

31.3,1.250000 

514 

7,1,1, -.000106 

570 


515 

8,2,1, -.000106 

571 

CONTINUE 

516 

9,3,1, -.000106 

572 

BOUNDARY CHANGE 

517 

THERMAL LOADS 

573 

3, 

518 

1, 

574 

7,1,1, -.000108 

519 

-.7,1.250000 

575 

8,2,1, -.000108 

520 


576 

9,3,1, -.000108 

521 

CONTINUE 

577 

THERMAL LOADS 

522 

BOUNDARY CHANGE 

578 

1, 

523 

3, 

579 

36.5,1.250000 

524 

7,1,1, .000063 

580 


525 

8,2,1, .000063 

581 

CONTINUE 

526 

9,3,1, .000063 

582 

BOUNDARY CHANGE 

527 

THERMAL LOADS 

583 

3, 

528 

1. 

584 

7,1,1, -.000126 

529 

6,1,1.250000 

585 

8,2,1, -.000126 

530 


586 

9,3,1, -.000126 
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587 

THERMAL LOADS 

641 

CONTINUE 

588 

1. 

642 

BOUNDARY CHANGE 

589 

41.1,1.250000 

643 

3. 

590 


644 

7,1,1, -.000180 

591 

CONTINUE 

645 

8,2.1, -.000180 

592 

BOUNDARY CHANGE 

646 

9,3,1, -.000180 

593 

3, 

647 

THERMAL LOADS 

594 

7,1,1, -.000142 

648 

1, 

595 

8,2,1, -.000142 

649 

52.0,1.250000 

596 

9,3,1, -.000142 

650 


597 

THERMAL LOADS 

651 

CONTINUE 

598 

1, 

652 

BOUNDARY CHANGE 

599 

45.0,1.250000 

653 

3, 

600 


654 

7,1,1, -.000179 

601 

CONTINUE 

655 

8,2,1, -.000179 

602 

BOUNDARY CHANGE 

656 

9,3,1, -.000179 

603 

3. 

657 

THERMAL LOADS 

604 

7,1,1, -.000155 

658 

1. 

605 

8,2,1, -.000155 

659 

50.7,1.250000 

606 

9,3,1, -.000155 

660 


607 

THERMAL LOADS 

661 

CONTINUE 

608 

1, 

662 

BOUNDARY CHANGE 

609 

48.1,1.250000 

663 

3, 

610 


664 

7,1,1, -.000175 

611 

CONTINUE 

665 

8,2,1, -.000175 

612 

BOUNDARY CHANGE 

666 

9,3,1, -.000175 

613 

3. 

667 

THERMAL LOADS 

614 

7,1,1, -.000166 

668 

1, 

615 

8,2,1, -.000166 

669 

48.6,1.250000 

616 

9,3,1, -.000166 

670 


617 

THERMAL LOADS 

671 

CONTINUE 

618 

1. 

672 

BOUNDARY CHANGE 

619 

50.4,1.250000 

673 

3, 

620 


674 

7,1,1, -.000243 

621 

CONTINUE 

675 

8,2,1, -.000243 

622 

BOUNDiVRY CHANGE 

676 

9,3,1, -.000243 

623 

3. 

677 

THERMAL LOADS 

624 

7,1,1, -.000174 

678 

1, 

625 

8,2,1, -.000174 

679 

45.7,1.250000 

626 

9,3,1, -.000174 

680 


627 

THERMAL LOADS 

681 

CONTINUE 

628 

1, 

682 

BOUNDARY CHANGE 

629 

51.8,1.250000 

683 

3. 

630 


684 

7,1,1, -.000034 

631 

CONTINUE 

685 

8,2,1, -.000034 

632 

BOUNDARY CHANGE 

686 

9,3,1, -.000034 

633 

3. 

687 

THERMAL LOADS 

634 

7,1,1, -.000179 

688 

1, 

635 

8,2,1, -.000179 

689 

41.9,1.250000 

636 

9,3,1, -.000179 

690 


637 

THERMAL LOADS 

691 

CONTINUE 

638 

1, 

692 

BOUNDARY CHANGE 

639 

52.3,1.250000 

693 

3, 

640 


694 

7,1,1, -.000031 
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695 

8,2,1, -.000031 

749 

7.5,1.250000 

696 

9,3,1, -.000031 

750 


697 

THERMAL LOADS 

751 

CONTINUE 

698 

1. 

752 

BOUNDARY CHANGE 

699 

37.5,1.250000 

753 

3, 

700 


754 

7,1,1, .000108 

701 

CONTINUE 

755 

8,2,1, .000108 

702 

BOUNDARY CHANGE 

756 

9,3,1, .000108 

703 

3, 

757 

THERMAL LOADS 

704 

7,1,1, -.000028 

758 

1, 

705 

8,2,1, -.000028 

759 

.7,1.250000 

706 

9,3,1, -.000028 

760 


707 

THERMAL LOADS 

761 

CONTINUE 

708 

1, 

762 

BOUNDARY CHANGE 

709 

32.4,1.250000 

763 

3, 

710 


764 

7,1,1, -.000063 

711 

CONTINUE 

765 

8,2,1, -.000063 

712 

BOUNDARY CHANGE 

766 

9,3,1, -.000063 

713 

3, 

767 

THERMAL LOADS 

714 

7,1,1, -.000024 

768 


715 

8,2,1, -.000024 

769 

-6.1,1.250000 

716 

9,3,1, -.000024 

770 


717 

THERMAL LOADS 

771 

CONTINUE 

718 

1. 

772 

BOUNDARY CHANGE 

719 

26.8, 1.250000 

773 

3, 

720 


774 

7,1,1, .000569 

721 

CONTINUE 

775 

8,2,1, .000569 

722 

BOUNDARY CHANGE 

776 

9,3,1, .000569 

723 

3, 

777 

THERMAL LOADS 

CM 

7,1,1, -.000020 

778 

1, 

725 

8,2,1, -.000020 

7 79 

-12.9,1.250000 

726 

9,3,1, -.000020 

780 


727 

THERMAL LOADS 

781 

CONTINUE 

728 

1, 

782 

BOUNDARY CHANGE 

729 

20.7,1.250000 

783 

3. 

730 


784 

7,1,1, .001229 

731 

CONTINUE 

785 

8,2,1, .001229 

732 

BOUNDARY CHANGE 

786 

9,3,1, .001229 

733 

3. 

787 

THERMAL LOADS 

734 

7,1,1, -.000015 

788 

1. 

735 

8,2,1, -.000015 

789 

-19.4,1.250000 

736 

9,3,1, -.000015 

790 


737 

THERMAL LOADS 

791 

CONTINUE 

738 

1. 

792 

BOUNDARY CHANGE 

739 

14.2,1.250000 

793 

3, 

740 


794 

7,1,1, .000086 

741 

CONTINUE 

795 

8,2,1, .000086 

742 

BOUNDARY CHANGE 

796 

9,3,1, .000086 

743 

3, 

797 

THERMAL LOADS 

744 

7,1,1, -.000011 

798 

1, 

745 

8,2,1, -.000011 

799 

-25.6,1.250000 

746 

9,3,1, -.000011 

800 


747 

THERMAL LOADS 

801 

CONTINUE 

748 

1. 

802 

BOUNDARY CHANGE 
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803 

3, 

857 

THERMAL LOADS 

804 

7,1,1, .000114 

858 

1, 

805 

8,2,1, .000114 

859 

-50.4,1.250000 

806 

9,3,1, .000114 

860 


807 

THERMAL LOADS 

861 

CONTINUE 

808 

1, 

862 

BOUNDARY CHANGE 

809 

-31.3,1.250000 

863 

3, 

810 


864 

7,1,1, .000224 

811 

CONTINUE 

865 

8,2,1, .000224 

812 

BOUNDARY CHANGE 

866 

9,3,1, .000224 

813 

3, 

867 

THERMAL LOADS 

814 

7,1,1, .000139 

868 

1. 

815 

8,2,1, .000139 

869 

-51.8,1.250000 

816 

9,3,1, .000139 

870 


817 

THERMAL LOADS 

871 

CONTINUE 

818 

1, 

872 

BOUNDARY CHANGE 

819 

-36.5,1.250000 

873 

3, 

820 


874 

7,1,1, .000230 

821 

CONTINUE 

873 

8,2,1, .000230 

822 

BOUNDARY CHANGE 

876 

9,3,1, .000230 

823 

3, 

877 

THERMAL LOADS 

824 

7,1,1, .000162 

878 

1, 

825 

8,2,1, .000162 

879 

-52.3,1.250000 

826 

9,3,1, .000162 

880 


827 

THERMAL LOADS 

881 

CONTINUE 

828 

1, 

882 

BOUNDARY CHANGE 

829 

-41.1,1.250000 

883 

3, 

830 


884 

7,1,1, .000247 

831 

CONTINUE 

885 

8,2,1, .000247 

832 

BOUNDARY CHANGE 

886 

9,3,1, .000247 

833 

3, 

887 

THERMAL LOADS 

834 

7,1,1, .000183 

888 

1, 

835 

8,2,1, .000183 

889 

-52.0,1.250000 

836 

9,3,1, .000183 

890 


837 

THERMAL LOADS 

891 

CONTINUE 

838 

1, 

892 

BOUNDARY CHANGE 

839 

-45.0,1,250000 

893 

3, 

840 


894 

7,1,1, .000048 

841 

CONTINUE 

895 

8,2,1, .000048 

842 

BOUNDARY CHANGE 

896 

9,3,1, .000048 

843 

3, 

897 

THERMAL LOADS 

844 

7,1,1, .000200 

898 

1, 

845 

8,2,1, .000200 

899 

-50,7,1.250000 

846 

9,3,1, .000200 

900 


847 

THERMAL LOADS 

901 

CONTINUE 

848 

1, 

902 

BOUNDARY CHANGE 

849 

-48.?, 1.250000 

903 

3, 

850 


904 

7,1,1, .000047 

851 

CONTINUE 

905 

8,2,1, .000047 

852 

BOUNDARY CHANGE 

906 

9,3,1, .000047 

853 

3, 

907 

THERMAI, LOADS 

854 

7,1,1, .000214 

908 

1, 

855 

8,2,1, .000214 

909 

-48.6,1.250000 

856 

9,3,1, .000214 

910 






911 

CONTINUE 

966 

9,3,1, .000025 

912 

BOUNDARY CHANGE 

967 

THERMAL LOADS 

913 

3, 

968 

1, 

914 

7,1,1, .000045 

969 

-20.7,1.250000 

915 

8,2,1, .000045 

970 


916 

9,3,1, .000045 

971 

CONTINUE 

917 

THERMAL LOADS 

972 

BOUNDARY CHANGE 

916 

1, 

973 

3, 

919 

-45.7,1.250000 

974 

7,1,1, .000019 

920 


975 

8,2,1, .000019 

921 

CONTINUE 

976 

9,3,1, .000019 

922 

BOUNDARY CHANGE 

977 

THERMAL LOADS 

923 

3. 

978 

1, 

924 

7,1,1, .000043 

979 

-14.2,1.250000 

925 

8,2,1, .000043 

980 


926 

9,3,1, .000043 

981 

CONTINUE 

927 

THERMAL LOADS 

982 

BOUNDARY CHANGE 

928 

1, 

983 

3, 

929 

-41.9,1.250000 

984 

7,1,1, .000013 

930 


985 

8,2,1, .000013 

931 

CONTINUE 

986 

9,3,1, .000013 

932 

BOUNDARY CHANGE 

987 

THERMAL LOADS 

933 

3. 

988 

1, 

934 

7,1,1, .000039 

989 

-7.5,1.250000 

935 

8,2,1, .000039 

990 


936 

9,3,1, .000039 

991 

CONTINUE 

937 

THERMAL LOADS 

992 

BOUNDARY CHANGE 

938 

1. 

993 

3, 

939 

-37.5,1.250000 

994 

7,1,1, -.000106 

940 


995 

8,2,1, -.000106 

941 

CONTINUE 

996 

9,3,1, -.000106 

942 

BOUNDARY CHANGE 

997 

THERMAL LOADS 

943 

3, 

998 

1, 

944 

7,1,1, .000035 

999 

-.7,1.250000 

945 

8,2,1, .000035 

1000 


946 

9,3,1, .000035 

1001 

CONTINUE 

947 

THERMAL LOADS 

1002 

BOUNDARY CHANGE 

948 

1, 

1003 

3, 

949 

-32.4,1.250000 

1004 

7,1,1, .000063 

950 


1005 

8,2,1, .000063 

951 

CONTINUE 

1006 

9,3,1, .000063 

952 

BOUNDARY CHANGE 

1007 

THERMAL LOADS 

953 

3v 

1008 

1, 

954 

7,1,1, .000030 

1009 

6.1,1.250000 

955 

8,2,1, .000030 

1010 


956 

9,3,1, .000030 

1011 

CONTINUE 

957 

THERMAL LOADS 

1012 

BOUNDARY CHANGE 

958 

1, 

1013 

3, 

959 

-26.8,1.250000 

1014 

7,1,1, -.000569 

960 


1015 

8,2,1, -.000569 

961 

CONTINUE 

1016 

9,' 1, -.000569 

962 

BOUNDARY CHANGE 

1017 

Tni. 'AL LOADS 

963 

3, 

1018 

1, 

964 

7,1,1, .000025 

1019 

12.5, '..250000 

965 

8,2,1, .000025 

1020 




1021 

CONTINUE 
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APPENDIX 8. DATA INPUT FOR OPEN NONSYMMETRIC TMF LOOP 


THERMOMECHANICAL LOOP NO. 3 WITH FUNCTIONAL THEORY 
15000 4 9 9 2 10 

2 


BOUNDARY CONDITIONS 
0 1 1 
BOUNDARY CONDITIONS 
1 , 

l,3,l,l,l.E-8 
THERMAL LOADS 

I, 

0.,l.E-8 


PROPERTY 

1 

26.0E6 

1 4 

POST 

2 16 
1 

11 

PRINT CHOICE 


O.E-5 
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51 

1 

104 

7,1,1, -.000344 

52 

CONTROL 

105 

8,2,1, -.000344 

53 

49,15,1,0 

106 

9,3,1, -.000344 

54 

.05 

107 

THERMAL TOADS 

55 

RESTART 

108 

1. 

56 

1.1,0, 8, 8 

109 

68.9,2.104167 

57 

INITIAL STATE 

110 


58 

1,1,1 

111 

CONTINUE 

59 

1,4, 1,4 

112 

BOUNDARY CHANGE 

60 

940. 

113 

3, 

61 

END OPTION 

114 

7,1,1, -.000378 

62 

BOUNDARY CHANGE 

115 

8,2,1, -.000378 

63 

3, 

116 

9,3,1, -.000378 

64 

7,1,1, .000000 

117 

THERMAL LOADS 

65 

8,2,1, .000000 

118 

1, 

66 

9,3,1, .000000 

119 

120 
121 

79.6,2.104167 

67 

68 

THERMAL LOADS 
1. 

.0,2.104167 

CONTINUE 

69 

122 

BOUNDARY CHANGE 

70 


123 

3, 

71 

CONTINUE 

124 

7,1,1, -.000411 

72 

BOUNDARY CHANGE 

125 

8,2,1, -.000411 

73 

3, 

126 

9,3,1, -.000411 

74 

7,1,1, -.000085 

127 

THERMAL LOADS 
1, 

75 

8,2,1, -.000085 

128 

76 

9,3,1, -.000085 

129 

130 

131 

86.5,2.104167 

77 

78 

THERMAL LOADS 
1. 

19.5,2.104167 

CONTINUE 

79 

132 

BOUNDARY CHANGE 

80 


133 

3, 

%J\J 

81 

CONTINUE 

134 

7,1,1, -.000424 

82 

BOUNDARY CHANGE 

135 

8,2,1, -.000424 

83 

3^ 

136 

9,3,1, -.000424 

84 

7,1,1, -.000166 

137 

THERMAL LOADS 

85 

8,2,1, -.000166 

138 

1 

86 

9,3,1, -.000166 

139 

89.2,2.104167 

87 

THERMAL LOADS 

140 


88 

1, 

141 

CONTIwUE 

89 

38.1,2.104167 

142 

BOUNDARY CHANGE 

90 


143 

3. 

91 

CONTINUE 

144 

7,1,1, -.000401 

92 

BOUNDARY CHANGE 

145 

8,2,1, -.000401 

93 

3, 

146 

9,3,1, -.000401 

94 

7,1,1, -.000274 

147 

THERMAL LOADS 
1. 

95 

8,2,1, -.000274 

148 

96 

9,3,1, -.000274 

149 

87.6,2.104167 

97 

THERMAL LOADS 

150 

CONTINUE 

98 

1, 

54.8,2.104167 

151 

99 

152 

BOUNDARY CHANGE 

100 


153 

3 1 

101 

CONTINUE 

154 

7,1,1, -.000374 

102 

BOUNDARY CHANGE 

155 

8,2,1, -.000374 

103 

3. 

156 

9,3,1, -.000374 

157 

THERMAL LOADS 
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158 

1, 

211 

159 

81.8,2.104167 

212 

160 


213 

161 

CONTINUE 

214 

162 

BOUNDARY CHANGE 

215 

163 

3. 

216 

164 

7,1,1, -.000246 
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983 

3, 

928 

1. 

984 

7,1,1, .000225 

929 

-72.0,2.104168 

985 

E,2,l, .000225 

930 


986 

';,3,1, .000225 

931 

CONTINUE 

987 

THERMAL LOADS 

932 

BOUNDARY CHANGE 

988 

1, 

933 

3, 

989 

-68.9,2.104168 

934 

7,1,1, .000462 

990 


935 

8,2,1, .000462 

991 

CONTINUE 

936 

9,3,1, .000462 

992 

BOUNDARY CHANGE 

937 

THERMAL LOADS 

993 

3, 

938 

1. 

994 

7,1,1, .000070 

939 

-81.8,2.104166 

995 

8,2,1, .000070 

940 


996 

9,3,1, .000070 

941 

CONTINUE 

997 

THERMAL LOADS 

942 

BOUNDARY CHANGE 

998 

1, 

943 

3. 

999 

-54.8,2.104166 

944 

7,1,1, .000366 

1000 


945 

8,2,1, .000366 

1001 

CONTINUE 

946 

9,3,1, .000366 

1002 

BOUNDARY CHANGE 

947 

THERMAL LOADS 

1003 

3, 

948 

1, 

1004 

7,1,1, .000048 

949 

-87.6,2.104166 

1005 

8,2,1, .000048 

950 


1006 

9,3,1, .000048 

951 

CONTINUE 

1007 

THERMAL LOADS 

952 

BOUNDARY CHANGE 

1008 

1, 

953 

3, 

1009 

-38.1,2.104166 

954 

7,1,1, .000373 

1010 


955 

8,2,1, .000373 

1011 

CONTINUE 
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